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What can emission lines tell us? 



Emission lines are observed almost everywhere in the Universe, from the 
Earth's atmosphere (see Wyse & Gilmore 1992 for a summary) to the 
most distant objects known (quasars and galaxies), on all scales, and 
at all wavelengths, from the radio domain (e.g. Lobanov 2005) to the 
gamma rays (e.g. Diehl et al. 2006). They provide very efficient tools 
to explore the Universe, measure the chemical composition of celestial 
bodies, and determine the physical conditions prevailing in the regions 
where they are emitted. 

The subject is extremely vast. In these lectures, we will restrict our- 
selves in wavelength, being mostly concerned with the optical domain, 
with some excursions to the infrared and ultraviolet domains and, occa- 
sionally, to the X-rays. 

We will mainly deal with the mechanisms of line production, and with 
the interpretation of line intensities in different astrophysical contexts. 
We will not discuss quasars and Seyfert galaxies, since those are the 
subject of the lectures by Brad Peterson, nor Lyman a galaxies, which 
are extensively covered by Daniel Schaerer in this book. However, we 
will discuss diagnostic diagrams used to distinguish active galaxies from 
other emission line galaxies, and will mention some topics linked with 
H Lya. Most of our examples will be taken from recent literature on 
planetary nebulae, Hn regions and emission line galaxies. Emission line 
stars are briefly described in the lectures by Stephen Eikenberry in the 
same volume, and a more detailed presentation is given in the book "The 
Astrophysics of Emission Line Stars" by Kogure & Leung (2007) . 

The vast subject of molecular emission lines has been left aside. The 
proceedings of the symposium "Astrochemistry: Recent Successes and 
Current Challenges" (Lis et al. 2006) give a fair introduction to this 
rapidly expanding field. 
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In the present text, we will not go into the question of dopplcr shifts 
or line profiles, which tell us about radial velocities, and thus about 
dynamics. This is of course a very important use of emission lines, which 
would deserve another series of lectures. For example, in such objects 
as planetary nebulae or supernova remnants, emission line profiles allow 
one to measure expansion velocities in these objects, and thus investigate 
their dynamics. The distribution of radial velocities of planetary nebulae 
in galactic haloes is a way to probe their kinematics and infer the dark 
matter content of galaxies (Romanowsky 2006). Rcdshift surveys which 
are used to map the 3D distribution of galaxies in the Universe strongly 
rely on emission lines (e.g. Lilly et al. 2007), which are the safest tool 
to measure redshifts. 

We will however mention the great opportunity offered by integral- 
field spectrosocopy at high spectral resolution, which provides line in- 
tensities and profiles at every location in a given field of view. With 
appropriate techniques, this allows one to recover the 3D geometry of a 
nebula. 

The purpose of these lectures is not to review all the literature on 
ionized nebulae, but rather to give clues for understanding the infor- 
mation given by emission lines, to provide some tools for interpreting 
one's own data, and to argue for the importance of physical arguments 
and common sense at each step of the interpretation process. Therefore, 
we will review methods rather than objects and papers. This series of 
lectures complements in some sense the lectures entitled "Abundance 
determinations in Hn regions and planetary nebulae" (Stasihska 2004), 
where the reader is referred to. In order to save space, the topics that 
have been extensively treated there will not be repeated, unless we wish 
to present a different approach or add important new material. 

In the following, we will assume that the reader is familiar with Sect. 
1 (Basic physics of photoionized nebulae), Sect. 2 (Basics of abundance 
determinations in ionized nebulae) and Sect. 3 (Main problems and 
uncertainties in abundance determinations) of the lectures by Stasinska 
(2004). We also recommend reading Ferland's outstanding (2003) review 
"Quantitative spectroscopy of photoionized clouds". Those wishing a 
more complete description of the main physical processes occuring in 
ionized nebulae should consult the text books "Physical Processes in the 
Interstellar Medium" by Spitzer (1978), "Physics of Thermal Gaseous 
Nebulae" by Aller (1984), "Astrophysics of the Diffuse Universe" by 
Dopita & Sutherland (2003), and "Astrophysics of Gaseous Nebulae and 
Active Galactic Nuclei" by Osterbrock & Ferland (2006). For a recent 
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update on X-ray astrophysics, a field which is rapidly developing, one 
may consult the AIP Conference Proceedings on "X-ray Diagnostics of 
Astrophysical Plasmas: Theory, Experiment, and Observation" (2005, 
Volume 774, edited by Randall K. Smith). 

1.1 Generalities 

This section is meant to complement the Section "Basic physics of pho- 
toionized nebulae" , which presented the main physical processes occur- 
ing in ionized nebulae. 

1.1.1 Line production mechanisms 

Emission lines arise in diffuse matter. They are produced whenever an 
excited atom (or ion) returns to lower-lying levels by emitting discrete 
photons. There are three main mechanisms which produce atoms (ions) 
in excited levels: recombination, collisional excitation, and photoexcita- 
tion. 

1.1.1.1 Recombination 
Roughly two thirds of the recombinations of an ion occur onto excited 
states from which deexcitation proceeds by cascades down to the ground 
state. The resulting emission lines are called recombination lines and 
are labelled with the name of the recombined ion, although their intensi- 
ties are proportional to the abundance of the recombining species. The 
most famous (and most commonly detected) ones are H I lines (from the 
Balmer, Paschen, etc ... series) which arise from recombination of H + 
ions; Hei lines (A5876 ...) which arise from recombination of He + ; and 
Hen lines (A4686 ...) which arise from recombination of He ++ ions. Re- 
combination lines from heavier elements are detected as well, (e.g. C n 
A4267, On A4651 ...) but they are weaker than recombination lines of 
hydrogen by several orders of magnitude, due to the much lower abun- 
dances of those elements. 

The energy eiji emitted per unit time in a line I due the recombination 
of the ion j of an element X 1 can be written as 

e«, - n e n(X/)e° 7 T-«, (1.1) 

where is a constant and the exponent a is of the order of 1. Thus, 
recombination line intensities increase with decreasing temperature, as 
might be expected. 
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1.1.1.2 Collisional excitation 

Collisions with thermal electrons lead to excitation onto levels that are 
low enough to be attained. Because the lowest lying level of hydrogen 
is at 10.2 eV, collisional excitation of hydrogen lines is effective only at 
electron temperatures T e larger than ~ 2 x 10 4 K. On the other hand, 
heavy elements such as nitrogen, oxygen, neon have low lying levels 
which correspond to fine structure splitting of the ground level. Those 
can be excited at any temperature that can be encountered in a nebula, 
giving rise to infrared lines. At "typical" nebular temperatures of 8000 - 
12000 K, levels with excitation energies of a few eV can also be excited, 
giving rise to optical lines. Slightly higher temperatures are needed to 
excite levels corresponding to ultraviolet lines. 

In the simple two-level approximation and when each excitation is 
followed by a radiative deexcitation, the energy e^; emitted per unit 
time in a line I due to collisional excitation of an ion j of an element X % 
can be written as 

em = n e n{Xi) qijl hui = &.6310-% e n(X 3 i )Q ijl /uj ijl T- - 5 e^/ kT ^hu l , 

(1.2) 

where fiyj is the collision strength, Wyj is the statistical weight of the 
upper level, and Xiji is the excitation energy 

Collisionally excited lines (CELs) are traditionally separated into for- 
bidden, semi-forbidden and permitted lines, according to the type of 
electronic transition involved. Observable forbidden lines have transi- 
tion probabilites of the order 10~ 2 s _1 (or less for infrared lines), semi- 
forbidden lines of the order of 10 2 s , and permitted lines of the order 
of 10 8 s . This means that the critical density (i.e. the density at which 
the collisional and radiative deexcitation rates are equal) of forbidden 
lines is much smaller than that of intercombination lines and of per- 
mitted lines. Table 2 of Rubin (1989) lists critical densities for optical 
and infrared lines, while Table 1 of Hamman ct al. (2001) gives critical 
densities for ultraviolet lines. 

Resonance lines are special cases of permitted lines: they are the 
longest wavelenght line arising from ground levels. Examples of reso- 
nance lines are H Lya, Civ A1550. 

1.1.1.3 Fluorescent excitation 

Permitted lines can also be produced by photoexcitation due to stellar 
light or to nebular recombination lines. The Bowen lines (Bowen 1934) 
are a particular case of fluorescence, where Olll is excited by the Hen 
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Lya line, and returns to the ground level by cascades giving rise in 
A3133, 3444, as well as to the line Oili A304 which in turn excites a 
term of Nlll. 

The interpretation of fluorescent lines is complex, and such lines are 
not often used for diagnostics of the nebulae or their ionizing radiation. 
On the other hand, it is important to know which lines can be affected 
by fluorescence, in order to avoid improper diagnostics assuming pure 
recombination (see Escalante & Morisset 2005). Order of magnitude 
estimates can be done with the simple approach of Grandi (1976). 

Quantitative analysis of fluorescence lines requires heavy modelling. 
It can be used to probe the Hen radiation field in nebulae (Kastner & 
Bhatia 1990). The fluorescent X-ray iron line has been used to probe 
accretion flows very close to massive black holes (Fabian et al. 2000). 

1.1.1-4 Some hints 

Each line can be produced by several processes, but usually there is one 
which dominates. There are, however, cases when secondary processes 
may not be ignored. For example, the contribution of collisional ex- 
citation to H Lya is far from negligible at temperatures of the order 
of 2 x 10 4 K. The contribution of recombination to the intensities of 
[Oil] AA7320, 7330 is quite important at low temperatures (below, say, 
5000 K) and becomes dominant at the lowest temperatures expected in 
Hn regions (see Stasihska 2005). 

In the appendix, we give tables of forbidden, semi-forbidden and res- 
onance lines for ions of C, N, O, Ne, S, CI, Ar that can be observed in 
H 11 regions and planetary nebulae. These are extracted from the Atomic 
line list maintained by Peter van Hoof, available at 
Ihttp: / /www. pa. uky.edu / "peter / atomicl This database contains identi- 
fication of about one million allowed, intercombination and forbidden 
atomic transitions with wavelengths in the range from 0.5A to 1000/^m. 

W. C. Martin and W. L. Wiese produced a very useful atomic physics 
"compendium of basic ideas, data, notations and formulae" available at: 
|http://www.physics.nist.gov/Pubs/AtSpec/index.html| 

1.1.1.5 Atomic data 

In the interpretation of emission lines, atomic data play a crucial role. 
Enormous progress has been done in atomic physics these recent years 
but not all relevant data are yet available or known with sufficient ac- 
curacy. The review by Kallman & Palmeri (2007) is the most recent 
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critical compilation of atomic data for emission line analysis and pho- 
toionization modelling of X-ray plasmas. A recent assessment of atomic 
data for planetary nebulae is given by Bautista (2006). 
Many atomic data bases are available on internet. 

- A compilation of data bases for atomic and plasma physics is found at: 
|http://plasma-gate.weizmann.ac.il/ DBfAPP. html| 

- Reference data at: 

|http: / /physics.nist.gov/ PhysRefData/~physica] 

- Ultraviolet and X ray radiation at: 

http: / /www. arcetri.astro.it / science / chianti/chianti.html 

- Atomic data for astrophysics (but only up to 2000) at: 
ht t p ■ JJwww . pa . uky.edu / ~ verner / atom. html| 

- Atomic data from the Opacity Project at: 
http: / /cdsweb.u-strasbg.fr/topbase/topbase.html 

- Atomic data from the IRON project at: 
http: / / cdsweb.u-strasbg.fr/tipbase /home. html 



1.1.2 The transfer of radiation 

This section is not to provide a detailed description of radiative transfer 
techniques, but simply to mention the problems and the reliability of 
the methods that are used. 

1.1.2.1 The transfer of Lyman continuum photons emitted by the 

ionizing source 

The photons emitted by a source of radiation experience geometrical 
dilution as they leave the source. They may be absorbed on their way by 
gas particles, predominantly by hydrogen and helium. The first photons 
to be absorbed are those which have energies slightly above the ionization 
threshold, due to the strong dependence of the photoionization cross 
section with frequency (roughly proportional to v~ 3 ). This gives rise to 
a "hardening" of the radiation field as one approaches the outer edge 
of an ionized nebula. Photons may also be absorbed (and scattered) by 
dust grains. 

1.1.2.2 The transfer of the ionizing photons produced by the nebula 

Recombination produces photons that can in turn ionize the nebular 
gas. These photons are emitted in all directions, so that their transfer is 
not simple. Authors of photoionized codes have developed several kinds 
of approximations to deal with this. 
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The "on-the-spot approximation" (or OTS) assumes that all the pho- 
tons recombining to the ground-state are reabsorbed immediately and 
at the locus of emission. This is approximately true far from the ionizing 
source, where the population of neutrals is sufficiently large to ensure 
immediate reabsorption. However, in the zones of high ionization (or 
high "excitation" as often improperly said) this is not true. Compu- 
tationally, the OTS assumption allows one to simply discount all the 
recombinations to the excited levels. This creates a spurious tempera- 
ture structure in the nebula, with the temperature being overestimated 
in the high excitation zone, due to the fact that the stellar ionizing 
radiation held is "harder" than the combined stellar + recombination 
radiation held. The effect is not negligible, about 1000 - 2000 K in 
nebulae of solar chemical composition. Because of this, for some kind 
of problems, it might be preferable to use a simple ID photoionization 
code with reasonable treatment of the diffuse radiation than a 3D code 
using the OTS approximation. 

The OTS approximation, however, has its utility in dynamical simula- 
tions incorporating radiation transfer, because it significantly decreases 
the computational time. 

Note that the OTS approximation is valid on a global scale. In the 
integrated volume of a Stromgren sphere, the total number of ionizing 
photons of the source is exactly balanced by the total number of re- 
combinations to excited levels. This is a useful property for analytic 
estimations, since it implies that the total Ha luminosity of a nebula 
which absorbs all the ionizing photons (and is devoid of dust) is simply 
proportional to Qh, the total number of ionizing hydrogen photon£p. 

In increasing order of accuracy (and complexity) , then comes the "out- 
ward only approximation", which was hrst proposed in 1967 by Tarter 
in his thesis, very early in the era of photoionization codes (see Tarter ct 
al. 1969 for a brief description). Here, the ionizing radiation produced in 
the nebula is computed at every step, but artificially concentrated in the 
outward-directed hemisphere, where it is distributed isotropically. This 
gives a relatively accurate description of the nebular ionizing radiation 
held, since photons that are emitted inwards tend to travel without 
being absorbed until they reach the symmetrical point relative to the 
central source. The great advantage of this approximation is that it 

f This property is formally true only for a pure hydrogen nebula, but it so happens 
that absorption by helium and subsequent recombination produces line photons 
that ionize hydrogen and compensate rather well for the photons absorbed by 
helium 
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allows the computation of a model without having to iterate over the 
entire volume of the nebula. The code PHOTO used by Stasinska, and 
the code NEBU used by Pequignot and by Morisset are based on this 
approximation. The code CLOUDY, by Ferland, uses the outward-only 
approximation in a radial-only mode. It appears, from comparisons of 
benchmark models (e.g. Pequignot et al. 2001), that the global results 
of models made with codes that treat the transfer of diffuse radiation 
completely (e.g. the code nebula by Rubin) are quite similar. Note that 
the full outward-only approximation allows one to compute the ionizing 
radiation field in the shadows from optically thick clumps by artificially 
suppressing the stellar radiation field blocked by the clump. 

Codes treating the transfer of ionizing continuum photons exactly, it- 
erating over the entire nebula, interestingly appeared also at the begin- 
ning of the era of photoionization codes (Harrington 1968, Rubin 1968). 
At that time, computers were slow and had little memory, and only 
spherical or plane parallel geometries could be treated by such codes. 

With the present computational capacities, one can do much better 
and treat the transfer problem accurately for any geometry, by using 
Monte-Carlo methods. The first such code is MOCASSIN, by Ercolano 
(see Ercolano et al. 2003). One advantage of Monte-Carlo methods is 
that they allow one to treat the transfer accurately also in extremely 
dusty nebulae (Ercolano et al. 2005), while the outward-only approxi- 
mation breaks down in such cases. 

The transfer of resonance line radiation produced in the nebula is the 
most difficult to treat accurately, at least in classical approaches of the 
transfer. This is because it is generally treated in the "escape proba- 
bility" approximation. The effect of line transfer is crucial in optically 
thick X-ray plasmas such as the central regions of active galactic nuclei 
(AGN). The code TITAN by Dumont treats the transfer of line radi- 
ation in an "exact" manner, using the "Accelerated Lambda Iteration" 
method (Dumont et al. 2003). 

1.1.2.3 The non-ionizing lines emitted by the nebula 

In general, non-ionizing photons in dust-free nebulae escape as soon as 
they have been emitted (except perhaps in AGNs were column densities 
are higher). Resonance lines constitute one exception: they may be 
trapped a long time in the nebula, due to multiple scattering by atoms 
which, in nebular conditions, are predominantly in their ground level. 

In dusty objects, line photons suffer absorption and scattering by dust 
grains, on their path out of the nebula. Resonance lines, whose path 
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length can be multiplied by enormous factors due to atomic scattering, 
are then preferentially affected by dust extinction. Their observed lumi- 
nosities then represent only a lower limit to the total energy produced 
by these lines. 

Another case where the diagnostic potential is expected to be reduced 
due to transfer effects is that of infrared fine structure lines of abundant 
ions, such as [Oin] A88, 52 /im, which can become optically thick in 
massive Hn regions (Rubin 1978). However, due to a combination of 
independent reasons, this appears not to be the case even in the extreme 
situations explored by Abel et al. (2003). 

1.2 Empirical diagnostics based on emission lines 

1.2.1 Electron temperature and density 

It is well known, and mentioned in all textbooks, that some line ratios 
(e.g. the ratios of the lines labeled Al and N3 in Table 1.11) are strongly 
dependent on the temperature, since they have different excitation ener- 
gies. If the critical densities for collisional deexcitations are larger than 
the density in the medium under study, these line ratios only depend 
on the temperature, and are ideal temperature indicators. The most 
frequently used is the [Oni] A4363/5007 ratio. 

On the other hand, collisionally excited lines that arise from levels 
of similar excitation energies have their ratios that depend only on the 
density. The most common density indicator in the optical is the [S n] 
A6716/6731 ratio. Other ones can be easily found by browsing Table 
1.11. Rubin (1989) gives a convenient list of optical and infrared line 
density indicators showing the density range where each of them is use- 
ful. 

Similar plasma diagnostics are now available in the X-rays (Porquet 
& Dubau 2000, Delahaye et al. 2006, see also Porter & Fcrland 2006). 

1.2.2 Ionic and elemental abundances 

There are basically four methods to derive the chemical composition 
of ionized nebulae. The first one, generally thought to be the royal 
way, is through tailored photoionization modelling. The second is by 
comparison of given objects with a grid of models. These two methods 
will be discussed in the next section. In this section, which deals with 
purely empirical methods, we will discuss the other two: direct methods, 
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which obtain an abundance using information directly from the spectra, 
and statistical methods, which use relations obtained from families of 
objects. 

1.2.2.1 Direct methods 

In such methods, one first derives ionic abundance ratios directly from 
observed line ratios of the relevant ions: 

hjl = J n(X 3 i )n c € ij i(T c ,n c )dV ^ 
h'j'V J n(Xl)n e eiijn>(T e ,n e )dV 

where the Ps are the intensities and the e's are given by ey; = 
eijin(Xf)n e . 
Therefore: 

jn{X()n c dV = Im/hrv (1 4 , 

]n{X?,)n e dV em(Juni)le lly v{Tv,n v y 

where 7] and n/ are respectively the electron temperature and density 
representative of the emission of the line I. 

Assuming that the chemical composition is uniform in the nebula, one 
obtains the element abundance ratios: 

n(X v ) ey7(T i ,n;)/e i ' J vi'(T//,n;/) 
where icf is the ionization correction factor, equal to: 

JnjXiynjX^dV 
$n{X{,)/n{X v )n c dV' 

In the case where several ions of a same element are observed, one can 
use a "global" icf adapted to the ions that are observed (e.g. ic/(0 + + 
++ ) for planetary nebulae in which oxygen may be found in higher 
ionization stages) . Note that in H II regions (except those ionized by hot 
Wolf-Rayet stars), icf(0++ 0++)=l. 

The application of direct methods requires a correct evaluation of the 
Xj's and n;'s as well as a good estimate of the ionization correction factor. 

For some ions, the 7}'s can be obtained from emission line ratios 
such as [Om] A4363/5007 or [Nil] A5755/6584. For the remaining ions, 
the T;'s are derived using empirical relations with T (4363/5007) or 
T(5755/6584) obtained from grids of photoionization models. The most 
popular empirical relations are those listed by Garnett (1992). A newer 
set of relations, based on a grid of models that reproduces the properties 



12 



What can emission lines tell us? 



of Hn galaxies, is given by Izotov et al. (2006). It must be noted, how- 
ever, that observations show larger dispersion about those relations than 
predicted by photoionization models. It is not clear whether this is due 
to underestimated observational error bars, or to additional processes 
not taken into account by photoionization models. At high metallicities 
0, the relevance of any empirical relation between the various Ti 's is even 
more questionable, due to the existence of large temperature gradients 
in the nebulae, which are strongly dependent on the physical conditions. 

The ionization correction factors that are used are either based on 
ionization potential considerations, or on formulae obtained from grids 
of photoionization models. For Hii galaxies, a set of icf's is given by 
Izotov et al. (2006). For planetary nebulae, a popular set of ionization 
correction factors is that from Kingsburgh & Barlow (1994), which is 
based on a handful of unpublished photoionization models. Stasihska 
(2007, in preparation) gives a set of icfs for planetary nebulae based on 
a full grid of photoionization models. It must be noted, however, that 
theoretical icf's depend on the model stellar atmospheres that are used 
in the photoionization models. Despite the tremendous progress in the 
modelling of stellar atmospheres in the recent years, it is not yet clear 
whether predicted spectral energy distributions (SEDs) in the Lyman 
continuum are correct. 

Finally, note that the line-of-sight ionization structure, in the case 
of observations that sample only a small fraction of the entire nebula, 
is different from the integrated ionization structure. This is especially 
important to keep in mind when dealing with trace ionization stages. 

A case of failure of T e -based abundances: metal rich giant H II 
regions Until recently, it was not possible to measure the electron 
temperature in metal-rich H n regions. The usual temperature diag- 
nostics involve weak auroral lines which easily fall below the detec- 
tion threshold at low temperatures. With very large telescopes such 

f Throughout the paper the word "metallicity" is used with the meaning of "oxygen 
abundance". This is common practise in nebular studies. Although oxygen is 
not a metal according to the definition given by chemistry, in nebular astronomy 
the word metal is often used to refer to any element with atomic weight > 12. 
The use of the O/H abundance ratio to represent the "metallicity" - as first 
done by Peimbert (1978) — can be justified by the fact that oxygen represents 
about half of the total mass of the "metals", and that it is the major actor - 
after hydrogen and helium - for the emission spectrum of nebulae. Note that, for 
stellar astronomers, the word "metallicity" is related to the iron abundance, and 
not to the oxygen abundance, so that the two uses of the word "metallicity" are 
not strictly compatible, since the ratio O/Fe changes during the course of chemical 
evolution. 
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temperature-sensitive line ratos as [Oin] A4363/5007, [Nil] A5755/6584, 
[S in] A6312/9532 can now be measured even at high metallicities (e.g. 
Kcnnicutt et al. 2003, Brcsolin et al. 2005, Brcsolin 2007). How- 
ever, due to the large temperature gradients expected to take place in 
high-metallicity nebulae, which are a consequence of the extremely ef- 
ficient cooling in the ++ zone due to the infrared [O in] lines, [O in] 
A4363/5007 does not represent the temperatures of the ++ zone. As a 
consequence, the derived abundances can be strongly biased, as shown 
by Stasihska (2005). The magnitude of the bias depends on the physical 
properties of the Hn region and on which observational temperature 
indicators are available. 

A further problem in the estimation of T e at high metallicity is the 
contribution of recombination to the intensities of collisionally excited 
lines, which becomes important at low values of T e . For example, the 
contribution of recombination from ++ to the intensity of [O n] can be 
very important. It can be corrected for using the formula given in Liu et 
al. (2000), provided that the temperature characteristic of the emission 
of the recombination line is known. If the temperature is measured using 
ratios of CELs only, such is not the case. 

1.2.2.2 Statistical methods 

In many cases, the weak [O in] A4363 or [N n] A5755 lines are not avail- 
able because either the temperature is too low, or the spectra are of low 
signal-to-noise, or else the data consist of narrow band images in the 
strongest lines only. Then, one may use the so-called "strong line meth- 
ods" to derive abundances. Such methods are only statistical, in the 
sense that they allow one to derive the metallicity of an H II region only 
on the assumption that this Hn region shares the same properties as 
those of the H II regions used to calibrate the method. In practise, such 
methods work rather well for giant H n regions, since it appears that 
giant Hii regions form a narrow sequence (see e.g. McCall et al. 1985), 
in which the hardness of the ionizing radiation field and the ionization 
parameter are closely linked to the metallicity. Indeed, an increased 
metallicity enhances the metal line blocking of the emergent stellar flux 
in the extreme ultraviolet and softens the ionizing spectrum. In addi- 
tion, the pressure exerted on the nebular gas increases with the strength 
of the stellar winds, which are related to metallicity, and this in turn 
decreases the ionization parameter (Dopita et al. 2006). 

Unlike direct methods for abundance determinations, statistical meth- 
ods have to be calibrated. The reliability of these methods depends not 
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only on the choice of an adequate indicator, but also on the quality of the 
calibration. This calibration can be done using grids of ab-initio pho- 
toionization models (McGaugh 1991), using a few tailored photoioniza- 
tion models (Pagel et al. 1979), abundances derived from direct methods 
(Pilyugin & Thuan 2005), or objects other than Hn regions thought to 
have the same chemical composition (Pilyugin 2003) 

The oldest and still most popular statistical method is the one based 
on oxygen lines. Pagel et al. (1979) introduced the ([O n] A3727 + [O III] 
A4959,5007)/H/3 ratio (later referred to as R23 or O23) to estimate O/H. 
This method has been calibrated many times, with results that may 
differ by about 0.5 dex. Mc Gaugh (1994) and later Pilyugin (2000, 
2001) refined the method to account for the ionization parameter. 

Many other metallicity indicators have been proposed: [O ill] A5007/[Nil] 
A6584 (0 3 N 2 ) by Alloin et al. (1979); [Nil] A6584/H/3 (N 2 ), by Storchi- 
Bergmann et al. (1994) ; ([Sin] A9069+[Sn] A6716,6731)/Ha (S 23 ) by 
Vflchez&Esteban(1996); [Nn] A6584/[0n] A3727) (N 2 2 ) by Dopita et 
al. (2000); [Ann] A7135/[Om] A5007) (Ar 3 3 ) and [Sm] A9069/[Om] 
A3869) (S 3 3 ) by Stasihska (2006); [Nem] A9069/[On] A3727) (Ne 3 2 ) 
Nagao et al. (2006). The metallicity indicators proposed until 2000 
have been intercompared by Perez-Montero & Diaz (2005). However, 
all those methods will have to be recalibrated when the emission line 
properties of the most metal rich Hn regions are well understood, which 
is not the case at present. 

A few comments are in order. First, any method based on the ra- 
tio of an optical CEL and a recombination line (e.g. 23 or S 23 ) is 
bound to be double- valued, as illustrated e.g. by Fig. 7 of Stasihska 
(2002). This is because, at low metallicities, such ratios increase with 
increasing metallicity but, at high metallicities, they decrease due to the 
increased cooling by infrared lines, which lowers the temperature below 
the excitation threshold of optical CELs. In such circumstances, exter- 
nal arguments must be found to hnd out whether the object under study 
is on the "high abundance" or "low abundance" branch. The most com- 
mon argument is based on the [N 11] A6584 line. The reason why this 
argument works is that N/O is observed to increase as O/H increases, 
at least at high metallicities. Besides, high metallicity H II regions tend 
to have lower ionization parameters, favouring low-excitation lines such 
as [Nil] A6584. The biggest problem is at intermediate metallicities, 
where the maximum of 23 and S 23 occurs, and the metallicity is very 
ill-determined. By using both 23 and S 23 indexes at the same time, it 
would perhaps be possible to reduce the uncertainty. 
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Methods that use the [Nil] A6584 lines have another potential diffi- 
culty. Chemical evolution of galaxies changes the N/O ratio in a compli- 
cated and non universal way. Therefore, a calibration is not necessarily 
relevant for the group of objects under study. 

Perhaps the most satisfactory methods, on the theoretical side, are 
the ones using the Ar3C>3 or S3O3 indicators, since these indicators 
are monotonic and work for well understood reasons, directly linked 
to metallicity. 

On the contrary, the Ne302 index, which is seen to decrease as metal- 
licity decreases, behaves in such a way only because more metal rich 
giant Hll regions happen to be excited by a softer radiation field and 
have a lower ionization parameter. This is a very indirect metallicity 
indicator! 

It is important to be aware that, in principle, strong line methods 
can be safely used only when applied to the same category of objects 
that were used for the calibration. The meaning of the results in the 
case of integrated spectra of galaxies, for example, is far from obvious 
in an absolute sense. Such spectra contain the light from Hll regions 
differing in chemical composition and extinction as well as the light from 
the diffuse ionized interstellar medium. In addition, inclination effects 
may be important. A few studies have adressed these issues from an 
observational point of view (Zaritsky et al. 1994, Kobulnicky et al. 1999, 
Moustakas & Kennicutt 2006), but clearly the subject is not closed. 

A further step in strong line abundance determinations has been made 
by using ratios of line equivalent widths instead of intensities (Kobul- 
nicky et al. 2003). The advantage of using equivalent widths is that 
they are almost insensitive to interstellar reddening, which allows one 
to apply the method even when reddening corrections are not available, 
especially at redshifts larger than 1.6. The reason why equivalent widths 
work well for integrated spectra of galaxies is that there is empirically 
a very close correlation between line intensities and equivalent widths, 
meaning that, statistically, stellar and nebular properties as well as the 
reddening are closely interrelated. 

1.2.3 Estimation of the effective temperature of the ionizing 

stars 

T± from the Zanstra method. This method, proposed by Zanstra 
(1931), makes use of the fact that the number of stellar quanta in the 
Lyman continuum, normalized to the stellar luminosity at a given wave- 
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X (A) X (A) 

Fig. 1.1. Spectral energy distributions (SEDs) for effective temperatures corre- 
sponding to massive stars (left) and central stars of planetary nebulae (right). 
The dotted line indicates the position of the ionization potential of hydrogen, 
the full line indicates the wavelength of the V filter. 



length, is an increasing function of the effective temperature. This is 
illustrated in Fig. 1.1 (based on modern stellar model atmospheres). 
In practise, it is the luminosity of the H/3 line which is the counter of 
Lyman continuum photons. This assumes that all the Lyman contin- 
uum photons are absorbed by hydrogen. This assumption breaks down 
in the case of density bounded nebulae or of nebulae containing dust 
mixed with the ionized gas. In real nebulae, some Lyman continuum 
photons are also absorbed by He and He + . However, recombination 
of these ions produce photons that are able to ionize hydrogen, so that 
the basic assumption of the Zanstra method is generally remarkably well 
fulfilled. Of course, the value of the effective temperature, T*, obtained 
by the Zanstra method will depend on the model atmosphere used in 
the derivation. 

For very hot stars, such as the central stars of planetary nebulae, one 
can also define a He + Zanstra temperature, using the Hen A4686 flux 
as a measure of the number of photons with energies above 54.4 eV. 

A few properties, illustrated by the models shown in Fig. 1.2, are 
worth of noting. 

• He I A5876/H/3 measures T* only in a small range of temperatures (T* 
< 40 kK for our models with a blackbody as seen in Fig. 1.2. a, this 
limit being of course different for realistic stellar atmospheres) . 

• Due to the competition between H° and He + to absorb photons with 
energies above 54.4 eV, Hen A4686/H/3 saturates at T± > 150 kK and 
depends on U at T* > 100 kK (Fig. 1.2.b). 

• He II A4686/H/3 is independent of He/H (Fig. 1.2.b), while Hell 
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Fig. 1.2. Sequences of photoionization models for spherical, uniform nebulae 
ionized by stars of different effective temperatures. Thin black curve: ionisa- 
tion parameter U=W~ 2 , He/H=0.1. Thicker red curve: U=10~ 3 , He/H=0.1. 
Thick green curve: (7=10 -2 , He/H=0.15. For simplicity, the models were run 
with a blackbody central star. 



A4686/Hei A5876 depends on He/H (Fig. 1.2.c). Indeed, Hen A4686 
is a counter of photons, as long as the He ++ Stromgren sphere is 
smaller than the H + Stromgren sphere. 
• The temperatures in the H° and He + zones are not equal (Fig. 1.2. 
f)- 



T* from the observed ionization structure. Since the ionization 
structure of a nebula depends on one can think of using the line ratios 
of two successive ions to infer (e.g. Kunze et al. 1996). However, 
in such an approach, the effect of the ionization parameter must be 
considered as well. 

To alleviate this problem, Vflchez & Pagel (1988) introduced of a 
"radiation softness parameter" , which they defined as: 

0+ /° ++ 

and expected to be independent of U in first approximation. This, how- 
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ever, was not sufficient for an accurate determination of (or even for 
ranking the effective temperatures of different objects). 

Morisset (2004) constructed a grid of photoionization models with 
various values of U (at various metallicities) to allow a proper estimate of 
T*. The grid was constructed using the WM-Basic model atmospheres 
of Pauldrach et al. (2001). 

T* from energy-balance methods. This method was first proposed 
by Stoy (1933). It makes use of the fact that the heating rate of a nebula 
is a function of the effective temperature. Since, in thermal equilibrium, 
the heating rate is compensated by the cooling rate, a measure of the 
cooling rate allows one to estimate Most of the cooling is done 
through collisionally excited lines. Therefore, an estimation of T* can 
be obtained from the formula: 

^LcEL/L(Ul3) = f{T*). (1-8) 

Pottasch & Preite-Martinez (1983) proposed a calibration of this method 
for planetary nebulae. 

A similar argument led Stasihska (1980) to propose a diagram to esti- 
mate the average effective temperature of the ionizing stars of giant H n 
regions. This diagram plots the value of [Oin] A4363/5007 as a func- 
tion of the metallicity O/H, for photoionization models corresponding 
to various values of T*. 

1.2-4 Determining the star formation rate 

The star formation rate is an important quantity to measure in galax- 
ies. It can be obtained in many ways: ultraviolet continuum, far infrared 
continuum, radio continuum, recombination lines, forbidden lines. The 
last two are sensitive to the most recent star formation (less that a few 
Myr, which is the lifetime of the most massive stars). Note that, while 
the luminosity in a line measures the absolute value of the recent star 
formation rate, the equivalent width of the line measures the ratio of the 
present to past star formation rate. Any technique must be calibrated 
using simulations of stellar populations in which the basic parameters 
are the stellar initial mass function, the star formation history, and the 
metallicity. These simulations are based on libraries of stellar evolution- 
ary tracks and of stellar atmospheres. Kennicutt (1998) and Schaerer 
(2000) give exhaustive reviews on the question. Here, we simply men- 
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[Oil] 3727 / [Oil] 3727 / [Oil] 3727 / H0 

Fig. 1.3. Diagrams using data on galaxies from the Sloan Digital Sky Survey 
(taken from Stasiriska 2006) to show how strongly the [O n] A3727/H/3 ratio 
varies among normal star forming galaxies (left), galaxies containing an active 
nucleus (right), and hybrid galaxies (middle). 

tion a few issues regarding the estimation of star formation rates using 
emission lines. 

Star formation rate using H/3 or Ha. As any determination of 
a physical parameter using H/3, the basic assumption is that all the 
Lyman continuum photons are absorbed by the gas. Also, the effect of 
intervening dust extinction must be properly accounted for (see Kewley 
et al. 2002). As shown by Schaerer (2000), the derived star formation 
rate strongly depends on the adopted stellar initial mass function and 
upper stellar mass limit. 

Star formation rate using [O n] A3727. It may seem strange to 
use such a line to measure the star formation rate, instead of simply 
using Ha. In addition to the provisos on Ha, the [Oil] A3727 obviously 
must depend on the metallicity and ionization parameter. The main 
reason for attempting to use the [O n] A3727 line in spite of this is that 
it can be observed in the optical range at larger redshifts than Ha or 
H/3. However, it is an observational fact that the [On] A3727/H/3 ratio 
strongly varies among emission line galaxies, even discarding objects 
containing an active nucleus (see Fig. 1.3). Therefore, the use of [On] 
A3727 as a star formation rate indicator is extremely risky. 

1.2.5 How to distinguish normal galaxies from A GN hosts ? 

After the discovery of spiral galaxies with a very bright nucleus emitting 
strong and broad (several thousands of km s _1 ) emission lines (Seyfert 
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1941), it became clear that these galactic nuclefj] were the locus of vi- 
olent, non-stellar activity (Burbidge et al. 1963, Osterbrock & Parker 
1965), perhaps of the same nature as found in quasars. Heckman (1980) 
performed a spectroscopic survey of the nuclei of a complete sample 
of 90 galaxies, and found that the presence of low ionization nuclear 
emission-line regions (LINERs) were quite common, and seemed to be 
the scaled down version of Seyfert nuclei. Baldwin et al. (1981) were the 
first to propose spectroscopic diagnostics based on emission line ratios 
to distinguish normal star forming galaxies from active galactic nuclei 
(AGN). The most famous is the [Om] A5007/H/3 vs [Nil] A6584/Ha 
diagram, often referred to as the BPT diagram (for Baldwin, Phillips, 
& Terlevich). The physics underlying such a diagram is that photons 
from AGNs are harder than those from the massive stars that power 
Hn regions. Therefore they induce more heating, implying that optical 
collisionally excited lines will be brighter with respect to recombination 
lines than in the case of ionization by massive stars only. Veilleux & 
Osterbrock (1987) proposed additional diagrams: [Om] A5007/H/3 vs 
[Sn] A6725/Ha, and [Om] A5007/H/3 vs [Oi] A6300/Ha. As previously 
found by McCall et al. (1985), giant Hll regions form a very narrow 
sequence in these diagrams. 

It was a great surprise, after the first thousands of galaxy spectra from 
the Sloan Digital Sky Survey (SDSS, York et al. 2000) were released, 
to find that a proper subtraction of the stellar continuum in galaxies 
(Kauffmann et al. 2003) allowed one to see a second sequence in the BPT 
diagram, in the direction opposite to that of the star-forming sequence. 
Thus, emission line galaxies in the BPT diagrams are distributed in two 
wings, which look like the wings of a flying seagull (see Fig. 1.4). 

Just a few years before, Kewley et al. (2001) had constructed a grid of 
photoionization models in order to determine a theoretical upper limit 
to the ionization by massive stars in the BPT diagram. This upper limit, 
later referred to as the "Kewley line" proved well to the right of the star 
forming wing from the SDSS. Kauffman et al. (2003) shifted this line to 
the left to define an empirical limit between normal star forming galaxies 
and AGN hosts (the "Kauffmann" line). Stasinska et al. (2006) found 
this limit to be still too "generous" , and proposed a more restrictive 
one, based on a grid of photoionization models aimed at reproducing 
the upper envelope of the left wing of the seagull. Its equation is: 



f The term Seyfert galaxy was used for the first time by de Vaucouleurs (1960). 
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[Nil] 6583 / Ha 

Fig. 1.4. Galaxies from the Sloan Digital Sky Survey in the BPT diagram are 
distributed in two wings: the wings of a flying seagull. 

y = (-30.787 + 1.1358a; + 0.27297a; 2 )tanli(5. 7409a;) - 31.093, (1.9) 

where y = log [Om] A5007/H/3 and x = log [Nil] A6584/Ha. 

Note that, among the three diagnostic diagrams proposed by Veilleux 
and Osterbrock (1989), the [O ill] A5007/H/3 vs [Nil] A6584/Ha is the 
most efficient, due to the increase of N/O in galaxies of high metallicity. 
Then, heating by an AGN boosts the [Nil] line and creates a clear sepa- 
ration of the two wings. One then understands why the [O in] A5007/H/3 
vs [On] A3727/H/3 diagram, used as a surrogate of the BPT diagram 
at redshifts of ~ 0.2 (Lamareille et al. 2004) is much less efficient in 
separating AGN hosts from normal star forming galaxies. 

Stasinska et al. (2006) also noted that a classification based on [N il] 
A6584/Ha only is also feasible, and has the advantage of being applicable 
to a much larger number of galaxies, since [O ill] A5007 and H/3 are not 
needed. 

Finally, these authors proposed another classification diagram, plot- 
ting D n (4000j3 as a function of the equivalent width of [O n] A3727 (the 
DEW diagram). The rationale for this diagram is that if in a galaxy 
with a large £> n (4000) emission lines can be seen, they must be due to 
another cause than photoionization by massive stars. The hope is to be 
able to use this diagram to distinguish AGNs at redshifts higher than 

t -Dn(4000) is the discontinuity observed at 4000A in the spectra of galaxies; it 
increases with stellar metallicity and age. 
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0.2. However, the fact that higher redshifts correspond to younger ages 
of the stellar populations is a issue that requires further investigation. 

To end up this section, let us remark that the term "LINER" is nowa- 
days often employed to designate those galaxies that are found to the 
lower right of the BPT diagram. This is an unfortunate deviation from 
the original meaning, since it is likely that, in those galaxies, the emis- 
sion lines do not come only from the nucleus, but from a much larger 
zone of the galaxy. The original "LINERS" discovered by Heckman 
would not appear in the BPT diagram due to the low luminosity of the 
active nucleus, unless there is significant line emission from the rest of 
the galaxy. 



1.3 Photoionization modelling 

We now turn to the "royal way" of analysing emission line spectra: 
photoionization modelling. We will show that photoionization modelling 
is an art which not only requires photoionization codes, but also a certain 
dose of common sense. 



1.3.1 Publicly available photoionization codes 

Most of the codes listed below have been intercompared at the Lexing- 
tion conference "Spectroscopic Challenges of Photoionized Plasmas" in 
2000. The results of the code comparisons can be found in Pequignot et 
al. 2001. 

1.3.1.1 ID photoionization codes 

• CLOUDY by Gary Ferland and associates, computes models for ionized 
nebulae and photo-dissociation regions (PDR) . It is regularly updated, 
well documented, and widely used. It is available at 

http:/ /www. nublado.org 

• MAPPINGS by Michael Dopita + Kewley, Evans, Groves, Sutherland, 
Bincttc, Allen, Leitherer computes models for photoionized nebulae 
and for planar shocks. It can be found at 

[http: / /www.ifa.hawaii.edu/~kewley/ Mappings 

• Xstar by Tim Kallman computes models for photoionized regions 
with special attention to the treatment of X-rays. It is found at 
http://heasarc.nasa.gov/lheasoft/xstar/xstar.htmll 
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1.3.1.2 3D photoionization codes 

• CLOUDY-3D by Christophe Morisset is a pseudo 3D code, based on 
CLOUDY. It allows quick modelling of 3D nebulae and visualization 
(including computation and visualization of line profiles). It can be 
found at 

|http://132.248.1.102/Cloudy^3D| 

• MOCASSIN by Barbara Ercolano is a full 3D Monte-Carlo photoioniza- 
tion code that also treat dust transfer in an accurate manner. It is 
available from 

be@star.ucl.ac.uk (see also |http://hea- www.harvard.edu/~bercolano l 

1.3.1.3 Other codes 

Many other, independent photoionization codes are mentioned in the 
literature (some of them benchmarked in Pequignot et al. 2001), but 
are not of public access. This is especially the case of hydrodynamical 
codes that include the physics of photoionization and the computation 
of line emission. In such codes, of course, the physics of radiation is 
treated in a simplified manner, since a simultaneous treatment of the 
macro- and microphysics requires tremendous computing power. 

1.3.2 Why do photoionization models? 

There can be many different reasons for which to build photoionization 
models. For example, one might want to: 

• Check the sensitivity of input parameters on observable properties. 

• Compute a grid of models for easy interpretation of a certain class of 
objects. 

• Calculate ionization correction factors. 

• Derive the chemical composition of a given nebula. 

• Estimate characteristics of the ionizing source. 

• Probe stellar atmosphere model predictions in the far ultraviolet. 

1.3.3 How to proceed? 

Each of the problems above requires a specific approach. It is always 
worth spending some time to find the best way to achieve one's goal. 
For example, if one wants to derive the chemical composition of a nebula 
by means of emission line fitting, it is not sufficient to find one solution. 
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One must explore the entire range of possible solutions, given the obser- 
vational constraints. This is not always easy. 

If the aim is to interpret one object, or a given class of objects, the first 
step is to collect all the observational constraints needed for this purpose. 
This includes monochromatic images as well as line intensities in various 
wavelengths ranges and in different apertures. Also it is important to 
characterize the ionizing sources as best as possible from the observations 
(visual magnitude, spectral type in the case of a single star, age of the 
ionizing stellar population in the case of a large collection of coeval 
stars). 

Then, one must define a strategy: How to explore a parameter space? 
How to deal with error bars? How to test the validity of a model? 

Finally, one must evaluate the result of the investigation: Was the goal 
achieved? If not, what does this imply? In fact, this last aspect is often 
overlooked, while it is potentially very instructive and is an incentive for 
progress in the field. 

In the following we present some commented examples. 

1.3-4 Abundance derivation by tailored model fitting 

The general procedure would be something like: 

1) Define the input parameters: The characteristics of the ionizing 
radiation field (luminosity, spectral energy distribution); The density 
distribution of the nebular gas; The chemical composition of the nebular 
gas; The distance. 

2) Use an appropriate photoionization code and compute a model. 

3) Compare the outputs with the observations (corrected for extinc- 
tion by intervening dust): The total observed Ha flux; The Ha surface 
brightness distribution and the angular size of Ha emitting zone; The 
line intensities etc... 

4) Go back to 1) and iterate until the observations are satisfacto- 
rily reproduced. Here "satisfactorily" means that all the observational 
data are reproduced within acceptable limits, those limits taking into 
account both the observational errors and the approximations of the 
model. These "limits" should be set by a critical analysis of the situa- 
tion, before running the models. It may be that one finds no solution. 
This is by no means a defeat. It is actually an important result too, 
which tells something about the physics of the object. But, in order to 
be useful, such a result must be convincing, in other words the reason 
why no solution can be found must be clearly explained. 
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For a good quality model fitting, it is important to: 

1) Use as many observational constraints as possible (not only line 
intensity ratios). 

2) Keep in mind that the importance of the constraint has nothing to 
do with the strengths of the lines. For example Hen A4686/H/3, which is 
of the order of a few percent at most in H n regions, indicates the pres- 
ence of photons above 54.4 eV, which are not expected in main sequence, 
massive stars (unless they are part of a X-ray binary system) ; [O in] 
A4363/5007, [Nm] A5755/6584 indicate whether the energy budget is 
well reproduced. It is too often ignored that photoionization models do 
not only compute the ionization structure of the nebulae, they also com- 
pute their temperature, which is essential in predicting the strengths of 
the emission lines. 

3) Recall that some constraints are not independent. For example, if 
[Om] A5007/H/? is fitted by the model, then [Om] A4959/H/3 should 
be fitted as well because the [Om] A5007/4959 ratio is fixed by atomic 
physics, since both lines originate from the same levejj]. In no case can 
the fact that both lines are fitted at the same time be taken as a success 
of the model. On the other hand, if they are not, this may indicate an 
observational problem, e.g. that the strong [O III] A5007 line is saturated. 
As a matter of fact, many observers precisely use this ratio as a check 
of the accuracy of their observations. 

4) Choose a good estimator for the "goodness of fit" , e.g. avoid using a 
X 2 minimization technique without being convinced that this is the most 
appropriate test in the case under study. All the observables should be 
fitted (within limits defined a priori). 

5) Try to visualize the comparison model-result as much as possible. 
Examples among many possibilities can be found in Stasinska & Schaerer 
(1999) or in Stasinska et al. (2004). 

Outcomes from model-fitting: A priori, the most satisfactory situ- 
ation is when all the observations are fitted within the error bars. This 
may imply (but not necessarily) that the model abundances are the real 
abundances. If the constraints are insufficient, the model abundances 
may actually be very different from the true ones. 

Quite often, some of the observations cannot be fitted. This means 
either that the observations are not as good as thought, or that the model 

f The lack of variation of the [Om] A5007/4959 ratio among quasars at various 
redshifts shows that the fine structure constant, a, does not depend on cosmic time 
(Bahcall et al. 2003). A far-reaching application of emission-line astrophysics! 
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does not represent the object well. Some assumptions in the modelling 
may be incorrect. For example the nebula has a different geometry than 
assumed, or the stellar ionizing radiation field is not well described, or 
an important heating mechanism is missing. In such a situation, the 
chemical composition is generally not known with the desired accuracy, 
a fact too often overlooked. 

An example of photoionization modelling without a satisfac- 
tory solution: The most metal-poor galaxy I ZW 18 (Stasihska & 
Schaerer 1999) 

The observational constraints were provided by narrow band imaging 
in nebular lines and stellar continuum, along with optical spectra giving 
nebular line intensities. The ionizing radiation field was given by a stellar 
population synthesis model aimed at reproducing the observed features 
from hot Wolf-Rayet stars. The conclusion from the modelling exercise 
was that, even taking into account strong deviations from the adopted 
spectral energy distribution of the ionizing radiation and the effect of 
postulated additional X-rays, the photoionization models yield too low 
a [Oni] A4363/5007 ratio by about 30%[j]. This significant discrepancy 
cannot be solved by expected inaccuracies in the atomic data. The 
missing energy is of the same order of magnitude as the one provided 
by the stellar photons. 

Interestingly, this paper was rejected by a first referee on the account 
that "the temperature of the nebula, as indicated by the ratio [O in] 
A4363/5007, is in all trial models out of the observed range by 30%. 
The authors seem to have tried everything possible [...] but still ... no 
secure results." A second referee accepted the paper without changes 
emphasizing that "The results are of highest interest for future stud- 
ies of ionized regions, mainly by the identification of a serious problem 
in the interpretation of the [Oiii] A4363/5007 ratio, which cannot be 
consistently reproduced by the models " ! 

Another example of photoionization modelling without a sat- 
isfactory solution: The giant H n region NGC 588 in the galaxy M33 
(Jamet et al. 2005). 

In this case, the observational constraints were provided by a detailed 

f Note that the observed geometrical constraint provided by the images was impor- 
tant in reaching this conclusion, since the authors also showed that assuming a 
simpler geometry with no central hole — but in complete disagreement with the 
observations - one could obtain a higher electron temperature, due to reduced H 
Lya cooling. 
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characterization of the ionizing stars (Jamet et al. 2004), narrow band 
imaging, long slit optical spectra, and far infrared line fluxes from the 
ISO satellite. In the analysis, aperture effects were properly taken into 
account. The temperature derived from the [Olll] A4363/5007 ratio is 
larger than the one derived from the [Oin] A88/i/5007 ratio by 3000K. 
It was found that no photoionization model could explain the observed 
temperature diagnostics. Unless the measured [O in] A88/1 flux is in error 
by a factor 2, this implies that the oxygen abundance in this object is 
uncertain by about a factor of two. 

Note that both this study and the former one used a ID photoioniza- 
tion code, while the study of the nebular surface brightness distribution 
indicates that the nebula is not spherical. The nebular geometry was 
taken into account in the discussion, and the conclusions reached in both 
cases are robust. 

An example of photoionization modelling with insufficient con- 
straints: The chemical composition of the Galactic bulge planetary 
nebula M 2-5. 

This is an object for which no direct temperature diagnostic was avail- 
able. Using photoionization modelling, Ratag (1992, thesis) and Ratag 
et al. (1997) claimed an oxygen abundance of one fourth solar. How- 
ever, as shown in Stasihska (2002), models fulfilling all the observational 
constraints can be constructed with oxygen abundances as different as 
0/H= I.2xI0" 3 and 2.4xI0~ 4 ! This is linked to the "double value 
problem" discussed in Sect 1.2.2.2. Because of only crude knowledge 
of the physics of nebulae and insufficient exploration of the parameter 
space, Ratag had reached a possibly erroneous conclusion. 

1.3.5 Abundance derivation using grids of models 

One is often interested in determining the metallicities of a large sam- 
ple of objects. The extreme case being emission line galaxies from the 
Sloan Digital Sky Survey, which are counted by tens of thousands. A 
tailored model fitting is out of reach in such cases. Apart from statistical 
methods, discussed in Sect. 1.2.2.2, one may consider building a vast, 
finely meshed grid of photoionization models and use them to derive the 
metallicities of observed emission line galaxies either by interpolation or 
by Bayesian methods. Such an approach has been used by Chariot & 
Longhetti (2000), or Brinchmann et al. (2004). This method is power- 
ful and appealing, since it can be completely automatized. However, it 
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is not the ultimate answer to the problem of abundance determination 
and it must be used wih some circumspection. As shown by Yin et al. 
(2007), when applied to a sample of galaxies in which the [Oin] A4363 
line is observed, allowing direct abundance determinations, this methods 
returns significantly higher mctallicitics than the direct method. Yin et 
al. argue that the reason lies in the abundances chosen in the model 
grid. Since the procedure uses simultaneously all the strong lines to de- 
rive O/H, any offset between the real N/O and the N/O adopted in the 
model grid induces an offset in the derived O/H. 

1.3.6 Testing model atmospheres of massive stars using Hu 

regions 

In an on-going study of the Galactic Hu region M43 and its ionizing 
star (Simon-Diaz et al., in preparation), one of the aims is to use the 
nebula to probe the spectral energy distribution modelled for its ioniz- 
ing star. The nebula has a rather simple structure for an Hu region, 
being apparently round. The characteristics of the ionizing star (effec- 
tive temperature and gravity) were obtained by fitting the stellar H and 
He optical lines using the stellar atmosphere code fastwind (Puis et 
al. 2005). The stellar luminosity was then obtained, knowing the dis- 
tance. The predicted stellar spectral energy distribution was used as an 
input to the photoionization code CLOUDY (great care was taken in the 
treatment of the absorption edges, which are not the same in CLOUDY 
and fastwind). The nebula was assumed spherically symmetric, the 
density distribution was chosen to reproduce the observed Ha surface 
brightness distribution. The starting nebular abundances were obtained 
using the classical T -based method. The observational constraints are 
the dereddened line intensities at various positions in the nebula. 

It was found that the distribution of emission line ratios across the 
nebula cannot be reproduced by the model. While this could imply 
that the tested model atmosphere does not predict the energy distribu- 
tion correctly, another hypothesis must be investigated first: that the 
geometry assumed for the nebula is wrong. In spite of its roundish ap- 
pearance, it is possible that M43 is not a sphere, but a blister seen face 
on. The constructed photoionization model would then correspond to its 
"spherical impostor", using the expression by Morisset et al. (2005). As 
shown by Morisset & Stasinska (2006), observation of the emission line 
profiles would allow one to distinguish between the two geometries. In 
the meantime, photoionization computations using the pseudo 3D code 
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CLOUDY_3D of Morisset (2006) is being carried out, to provide models 
for the blister geometry. 

1.3.7 Photoionization study of an aspherical nebula using a 
3D code: the planetary nebula NGC 1009 

The use of a 3D photoionization modelling is generally difficult, due to 
the number of parameters involved to describe the geometrical structure 
of a nebula and to the long computational time required for a model. 
A ID code is often sufficient to identify a physical problem that needs 
solution, as seen in Sect. 1.3.4. On the other hand, a 3D code is essen- 
tial when tackling a problem where the effect of geometry is important. 
Such is the case, for example, of the low-ionization knots observed in 
planetary nebulae, for which it has been said that they show nitro- 
gen enhancement by factors 2-5 (e.g. Balick et al. 1994). Subsequent 
studies argued that this conclusion might be wrong, due to the use of 
an improper ionization correction scheme (it was assumed that N/O = 
N + /0 + ). Goncalves et al. (2006) examined this problem in the case of 
NGC 7009, a planetary nebula with two conspicuous symmetrical knots. 
Their aim was to explore the possibility that the enhanced [Nil] emis- 
sion observed in the knots could be due to ionization effects, by building 
a 3D photoionization model of homogeneous chemical composition re- 
producing the observed geometry and spectroscopic pecularities. They 
simulated observations in different regions of the nebula, and found that 
the N+/0+ ratio varies strongly with position and can explain (at least 
partly) the apparent nitrogen enhancement in the knots. 

1.3.8 The interpretation of data from integral field 

spectroscopy 

Integral field spectroscopy is becoming a major tool to study nebulae 
and galaxies. A wealth of spatially resolved data will routinely become 
available. The question now is which tools to use to interpret such 
data. Obviously, starting by building 3D photoionization models is not 
a good way to proceed. As argued above, 3D modelling is difficult and 
time consuming, and a proper evaluation of the benefits of this kind 
of approach is recommended before such an endeavour. On the other 
side, one may be tempted to use results from grids of published ID 
photoionizaion models or observational diagnostic diagrams relevant for 
the integrated nebular light. This is risky. For example, one might 
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attribute the increase of [Nil] /Ha or [Sn]/Ha in certain zones to non- 
stellar ionization while it might merely be a linc-of-sight effect in a region 
close to an ionization front. 

1.4 Pending questions 

There are a few important aspects which have not been mentioned above, 
for the sake of a more linear presentation. Most of them have been 
extensively treated in Stasinska (2004), and will be mentioned here only 
briefly, with some updates if necessary. 

1.4-1 Correction for reddening, underlying stellar absorption 

and aperture effects 

Before being analyzed in terms of abundances or star formation rates or 
being compared to the results of photoionization models, the intensities 
of observed lines must be corrected for various effects. The presence 
of dust between the zone of emission and the observer attenuates the 
collected radiation and modifies its colour. It appears that the dust 
extinction curve, which was once considered universal, is actually a one- 
parameter function characterized by the value of the total-to-selective 
extinction R v = A V /E(B-V) (Fitzpatrick 1999, 2004). The canonical 
value of generally used for extinction corrections, is 3.1 or 3.2. 
However, the measurement of Rv using 258 Galactic O stars yields a 
distribution around this value with a dispersion of ± 0.5. Values of Rv 
as small as 1.6 or as large as 5 are found (Patriarchi et al. 2003). 

The extinction curve of a single star by interstellar dust is different 
from the obscuration curve of an extended nebula (Calzetti 2001), be- 
cause in the latter case the observed radiation includes some scattered 
light. From a sample of starburst galaxies observed from the far ultravi- 
olet to the near infrared, an obscuration curve for galaxies was derived 
(Calzetti 1997, Calzetti et al. 2000). 

The amount of reddening is usually determined by fitting the observed 
Balmcr decrement to the theoretical one at the appropriate temperature. 
In emission line galaxies, in giant Hil regions (and also in unresolved 
planetary nebulae), the observed Balmer lines are affected by the un- 
derlying assumption from the stellar component. This absorption can 
be determined empirically together with the extinction by fitting several 
Balmer lines to the theoretical Balmer decrement (Izotov et al. 1994). 
It can also be modelled, using spectral synthesis techniques to reproduce 



1.4 Pending questions 31 




D n (4000) D n (4000) 



Fig. 1.5. The stellar absorption equivalent withs at Ha (left) and H/3 (right) 
as a function of the discontinuity at 4000 A. They have been determined in a 
sample of galaxies extracted from the Digital Sky Survey by stellar population 
analysis using the method of Cid-Fernandes et al. (2005). The red curves 
correspond to the fits given by Eqs. 1.10 and 1.11. 



the stellar absorption features other than Balmcr lines, and reading out 
from the model spectrum that best fits the data what is the stellar 
absorption at the the Balmer lines wavelengths (Cid Fernandes et al. 
2005). This procedure applied to a sample of 20000 galaxies from the 
Sloan Digital Sky Survey (see Stasihska et al. 2006 for a description of 
the sample) shows a spectacular correlation between absorption equiv- 
alent widths at the Balmer line wavelengths and the discontinuity at 
4000A, 7J„(4000), as shown for the first time here (Fig. 1.5). Since 
-D n (4000) is relatively easy to measure, this provides an empirical way 
to estimate the underlying stellar absorption at Ha and H/3. The data 
shown in Fig. 5 can be modelled by the expressions: 

£W ab8 (Ha) = 9exp(-1.2z) + 1.3 (1.10) 

and 

£W abB (H/3) = 400exp(-4.x + 0.5) + 2.5, (1.11) 

where x = L>„(4000). 

Another issue is when the objects under study are more extended 
than the observing beam. This fact must be taken into account in the 
analysis. For example, spectra obtained with several instruments must 
be carefully combined. This is especially important - and difficult - 
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when merging ultraviolet or infrared data with optical ones, since in 
that case flux calibration problems may arise. In case of observations 
where a large number of lines is available, it is theoretically possible to 
combine all the data by comparing each observation to the results of 
a tailored photoionization model seen through an appropriate aperture. 
Most of the time, however, aperture corrections are performed in an 
approximate way, neglecting the effect of the ionization structure of the 
nebulae on the line fluxes recorded through small apertures. 

A further problem, when dealing with spectra of entire galaxies is the 
physical meaning of the derived abundance, especially in the presence 
of abundance gradients. This issue, which is also relevant to the simple 
"strong line indicators" approach, has not yet been fully addressed. 

1.4-2 Escape of ionizing radiation 

Most nebular studies assume that the nebulae are ionization bounded. 
This affects such issues as the estimate of stellar temperatures via the 
Zanstra method, the estimation of star formation rates in galaxies, or 
the outcome of tailored photoionization modelling. 

There is, however, growing evidence that many nebulae are density 
bounded in at least some directions. This was known for a long time for 
planetary nebulae. It is now recognized to be the case also for many giant 
Hll regions (Beckman et al. 2002, Stasinska & Izotov 2003, Castellanos 
et al. 2003). 

1-4-3 The importance of stellar energy distribution 

The outcome of photoionization models depends critically on the adopted 
ionizing SED, especially as regards the nebular ionization structure. 
There are two main issues in this respect. One is the description of the 
ionizing spectral energy distribution of an individual star. The other is 
how one deals with the case of ionization by a group of stars. 

Model atmospheres for hot stars. The atmospheres of stars are gen- 
erally extended and experience important deviations from LTE. Hun- 
dreds of thousands of metallic lines produce line blanketing and line 
blocking effects, and are at the origin of radiative driven winds. Modern 
stellar atmosphere codes are able to handle these aspects, although each 
one works with some degree of approximation at least in one aspect. 
The most broadely used in photoionization modelling of H n regions 
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arc tlusty (Hubeny & Lanz 1995), WM-BASIC (Pauldrach et al. (2001), 
and CMFGEN (Hillier & Miller 1998). What matters for the modelling 
of ionized nebulae is the SED. There are sizeable differences among the 
predictions from these codes in the ionizing continuum. It is not clear, 
at present, which model predictions are best suited for photoionization 
modelling. A few studies (e.g. Morisset et al. 2004a and 2004b) have 
used observed Hll region to test the predictions of the model atmo- 
spheres. However, this is a difficult task, and the answer is not clear yet. 
Simon-Diaz et al. (2007) are attacking this problem anew, with both a 
theoretical and an observational approach. 

The central stars of planetary nebulae can reach much higher effec- 
tive temperatures than massive stars. Temperatures of over 2 x 10 5 /,K 
have been found. When they are close to becoming white dwarfs, these 
stars have tiny atmospheres for which the plane-parallel approximation 
is reasonable. A grid of plane-parallel line blanketed NLTE atmospheres 
has been computed by Rauch (2003) using the code pro2 (Werner & 
Dreizler 1993). At earlier stages of the planetary nebulae, the effect of 
winds are expected to be significant. They are important for planetary 
nebulae central stars of [WR] type, which share the same atmospheric 
properties as massive Wolf-Rayet stars. The same kind of approach as 
for H II regions has been used for planetary nebulae to test the ionizing 
SEDs predicted by model atmospheres (e.g. Stasihska et al. 2004 who 
used the Potsdam code for expanding atmospheres in non-LTE described 
in Hamann & Koesterke 1998). 

Some grids of model atmospheres computed with the above codes are 
available, for example at the CLOUDY website. 

Ionization by a group of stars. In the case of giant Hll regions, 
where the ionization is provided by many stars, one generally relies on 
SEDs predicted by stellar population synthesis codes. Starburst 99 (Lei- 
therer 1999, Smith et al. 2002) or pegase (Fioc & Rocca-Volmerange 
1997, 1999) are public-access codes able to do this. The integrated SED 
of an entire stellar population is obviously quite different from the SED 
of the component stars, but the quality of the SED for photoioniza- 
tion studies is related to the reliability of the predictions in the Lyman 
continuum of the stellar atmosphere codes used to compute the stellar 
libraries. 

There is, however, another important issue, pointed out by Cerviho 
et al. (2000). It is the fact that in real Hll regions, statistical sampling 
strongly affects the ionizing SED, and SEDs predicted by traditional 
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stellar population synthesis models may be far from real. Practical con- 
siderations as regards the outcome from photoionization models can be 
found in Cerviho et al. (2003) and a detailed discussion on uncertainties 
in stellar synthesis models and "survival strategies" is given by Cerviiio 
& Luridiana (2005). 

One of the yet unsolved problems is the question of the often observed 
He II emission in Hn galaxies. Traditional models (Stasihska & Izotov 
2003) have difficulties in accounting for it properly While there are 
several possible explanations (X-ray binaries, X-rays from hot star at- 
mospheres, etc.), sampling effects may also be an issue that needs to 
be considered. 

Finally, in the case of a group of ionizing stars, one might wonder 
if the effect on the Hn region strongly depends on whether the stars 
are concentrated or whether they are distributed in the nebula. This 
question is being adressed by Ercolano et al. (2007) 

1.4.4 Dust 

The question of the effect of dust in ionized nebulae has been amply 
reviewed in Stasihska (2004). Among publicly available photoionization 
codes, CLOUDY, mappings and MOCASSIN include a treatment of dust 
grains coupled with the gas particles. 

When it comes to entire galaxies, a proper treatment of dust should 
take into account the large scale geometrical distribution of dust inside 
the galaxies. Chariot & Longhetti (2001) did a first step by including 
a different prescription for the absorption of photons from H 11 regions 
and from older stars, as proposed by Chariot & Fall (2000). Panuzzo 
et al. (2003) implemented results from photoionization simulations with 
CLOUDY into models of dusty galaxies computed with the code Grasil2 
(Silva et al. 1998), which include a sound treatment of all the aspects of 
dust reprocessing, providing a more realistic description of the integrated 
spectrum of a galaxy. 

1-4-5 Temperature fluctuations and the ORL/CEL 

discrepancy 

Since the seminal work by Peimbert (1967), there exists the suspicion 
that elemental abundances derived from optical collisional emission lines 
(CELs) might be plagued by the presence of temperature fluctuations 
in nebulae. If such temperature fluctuations exist, then the highest 
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temperature zones will favour the lines of higher excitation, and as a 
result the temperature derived from [O in] A4363/5007 will overestimate 
the temperature of the ++ zone, and result in an understimate of 
elemental abundances with respect to hydrogen. 

Fourty years after, nebular specialists (see e.g. the recent reviews by 
Esteban 2002, Liu 2002, Fcrland 2003, Pcimbcrt & Peimbert 2003) are 
still debating about this issue, which is of fundamental importance for 
the reliability of abundance determinations of nebulae. We are talking 
here of an effect of typically 0.2 dex in metal abundances relative to 
hydrogen. One can summarize the debate in four questions: 

1) Is there evidence for the existence of temperature fluctuations at a 
level that significantly biases the abundances from CELs? 

2) If yes, is there a reliable procedure to correct for this bias? 

3) What would be the possible causes of such temperature fluctua- 
tions? 

4) Is the problem the same in H II regions and in planetary nebulae? 
At present, there are no generally accepted answers to these questions. 

In the following, we simply collect some elements for the discussion, by 
supplementing what can be found in the above reviews and in Stasihska 
(2004). 

1) Direct evidence for temperature fluctuations high enough to be of 
sizeable consequence for abundance determinations is scarce. It is pro- 
vided by Hubble Space Telescope temperature mapping of the Orion 
nebula (O'Dell et al. 2003) or spectroscopy (Rubin et al. 2003). Indi- 
rect evidence comes from the comparison of temperatures measured by 
various means (such as [Oin] A4363/5007 and the Balmer or Paschen 
jump). There are many examples, both in Hn regions and planetary 
nebulae, which point at a typical r.m.s. temperature fluctuation, t 2 , of 
0.04 (see references in the review quoted above). However, recent deter- 
mination of Balmer jump temperatures in H n galaxies (Guseva et al. 
2007) give t 2 close to zero, on average. Note that these latter measure- 
ments concern metal poor H n regions. Another piece of evidence comes 
from the discrepancy of abundances derived from ORLs and CELs for 
the same ion. In this case, the implied t 2 vary considerably, from less 
that 0.02 to almost 0.1, with the values for Hn regions being small, and 
those for planetary nebulae spanning a wide range. 

2) Pcimbcrt (1967) and Peimbert & Costero (1969) devised a scheme, 
based on Taylor series expansion of the analytic expressions giving the 
emissivities, to correct for the abundance bias, with the help of several 
simplifying assumptions. It can, however, be shown, on the basis of a 
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two-zone toy model (Stasiriska 2002), that a proper description of tem- 
perature fluctuations and any correction scheme for abundance determi- 
nations require at least three parameters, not just two as in Peimbert & 
Costero (1969). Besides, due to the restricted number of observational 
constraints, the application of the Peimbert & Costero scheme implies 
some assumptions which are incorrect, like the equality of t 2 (0 ++ ) and 
i 2 (H+). 

3) Some possible causes for the presence of temperature fluctuations 
have been listed by Torres-Pcimbcrt & Peimbert (2003) and Peimbert 
& Peimbert (2006). These include shocks, magnetic reconnection, dust 
heating, etc... However, the energy requirements to obtain t 2 values of 
about 0.04 are very large, typically a factor 2 with respect to the en- 
ergy gains from gas photoionization. Note that it is not so difficult to 
obtain a t 2 of 0.02-0.04 for the entire nebular volume, since the main 
cooling lines change across the nebulae. What is difficult is to obtain 
such a large t 2 in the ++ or + zones only, which is what would be 
implied by the ORL/CEL discrepancy. Using a "hot spot representa- 
tion" of temperature fluctuations, Binette et al. (2001) propose a sound 
quantification of energy requirements for the observed t 2 and find that 
"the combined mechanical luminosity of stellar winds, champagne flows 
and photoevaporation flows from proplyds is insufficient to account for 
the temperature fluctuations of typical Hn regions." Note that, at high 
metallicities, the natural temperature gradients in H n regions easily lead 
to a formal value of t 2 (0 ++ ) of 0.04 or even more, because of the effect 
of strong [Olll] A88, 52/im cooling. The difficulty is to obtain such a 
large t 2 (0 ++ ) at metallicities smaller than 12+log (O/H) = 8.7. 

4) There is now a large data base of high quality observations of 
metal recombination lines (see references in Liu 2006 for planetary neb- 
ulae and in Garcia- Rojas & Esteban 2006 for Hn regions). Analyzing 
the behaviour of the abundance discrepancy factors in both types of ob- 
jects Garcia-Rojas & Esteban (2006) conclude that the source for the 
ORL/CEL abundance discrepancy must be different in Hn regions and 
in planetary nebulae. Since the paper by Liu et al. (2000), it has been 
argued that, in planetary nebulae, the ORL/CEL abundance discrep- 
ancy is likely due to the presence of hydrogen-poor inclusions in the 
body of the nebula (see e.g. Liu 2002, 2003, 2006, Tsamis et al. 2004). 
These metal-rich clumps would be too cool to excite the CELs. Quite 
a small mass of matter in the clumps is sufficient to reproduce the ob- 
served ORL/CEL discrepancies. The nature of these clumps could be 
photoevaporating planetesimals. Tsamis et al. (2003) and Tsamis & 
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Pcquignot (2005) suggest that the ORL/CEL discrepancy in Hn re- 
gions might also be due to hydrogen poor inclusions (obviously of differ- 
ent astrophysical origin) . This point of view is not universally accepted 
(see Garcia- Rojas & Esteban 2006 for a different opinion). Yet, it seems 
that the galaxy chemical enrichment scenario proposed by Tenorio-Tagle 
(1996), in which the oxygen-rich products from supernova explosions end 
up in metal rich droplets falling back onto the galaxies after an excursion 
in the intergalactic medium, would be able to explain quantitatively the 
ORL/CEL abundance discrepancies observed in Hn regions (Stasihska 
et al. 2007). 

The debate on temperature fluctuations and ORL/CEL discrepancies 
is far from closed. Yet its outcome will have important implications on 
the use of emission lines in the derivation of elemental abundances in 
galaxies. Peimbert et al. (2006) are already proposing to use recombina- 
tion line abundances to calibrate the strong line methods for abundance 
derivations. They argue that recombination line abundances are not 
affected by any temperature bias, which is true in the case of a homo- 
geneous chemical composition of the nebulae. If abundance inhomo- 
geneities are present, one must investigate the astrophysical significance 
of the abundances measured. Stasihska et al. (2007) have examined this 
issue in the case of Hll regions. 

1.4-6 Shocks and related issues 

Photoionization is not the only process leading to the formation of emis- 
sion lines. Cooling flows (Cox & Smith 1976, Fabian & Nulsen 1977) 
and interstellar shocks can also produce them. Shocks are ubiquitous 
in galaxies, caused by jets (jets associated to brown dwarfs, protostars 
or massive black holes), winds (stellar winds, protostars, active galactic 
nuclei) or supersonic turbulence. Dopita & Sutherland (1995, 1996) pro- 
vide a grid of models of line emission produced by pure shocks. However, 
the efficiency of ionization by stellar photons is such that, whenever hot 
stars are present, they are likely to dominate the ionization budget of the 
surrounding gas. The same can probably be said for radiation arising 
from accretion onto a massive black hole. However, it is true that, in 
order to explain observed emission line ratios, the possible contribution 
of shock-heating to photoionization must be examined. 

The effects of shocks that leave a signature on emission line spectra 
in photoionized nebulae can be summarized as follows: 

1) Shocks generate compression and thus high gas densities. 
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2) This compression locally reduces the ionization parameter, and thus 
enhances the low ionization emission lines such as [O n] , [N n] , [S n] , and 
[Oi]. 

3) Shock heating produces very high temperatures (of the order of 
millions of K), which lead to collisional ionization, and the production 
of such ions as 3+ or He ++ . 

4) This high temperature gas emits X-rays, which contribute to the 
ionization in front of the shock and behind it. 

5) In the region between the X-ray emitting zone and the low ex- 
citation zone, the temperature of the gas passes through intermediate 
temperatures of several tens of thousands of K, boosting the emission of 
ultraviolet CELs and of auroral lines. The [Oin] A4363/5007 ratio can 
be significantly enhanced. 

However, it must be realized that effects 2, 3 and possibly 5 are spec- 
troscopically indistinguishable from the effects of photoionization by an 
external X-ray source. As for the enhancement of the low ionization 
lines, it can simply be due to a small ionization parameter without the 
necessity of invoking shock heating (see e.g. the explanation proposed 
by Stasihska & Izotov 2003 for the behaviour of the [Oi] lines in Hll 
galaxies) . 

Finally, dynamical heating and cooling as well as non equilibrium 
photoionization may be an issue in certain cases, when the dynamical 
times are shorter than the radiative time. Only a few works have started 
exploring these avenues, which require important computing power and 
therefore cannot use the same approach to real nebulae as the ones 
discussed in the present text. An enlightening introduction to the role 
of dynamics in photoionizcd nebulae is given by Hcnncy (2006) . 



1.5 Appendix: Lists of useful lines and how to deal with them 

In this appendix, we list all the forbidden, semi-forbidden and resonance 
lines from C, N, O, Ne, S, CI and Ar ions which have ionization poten- 
tials smaller than 126 eV. These lines as well as all the information on 
them were extracted from the Atomic Line List compiled by Peter van 
Hoof and accessible at http://www.pa.uky.edu/~peter/atomic Permit- 
ted lines are far too numerous to be listed there. They are available 
in the Atomic Line List. However it is more convenient, for a first ap- 
proach, to use the identifications made by Zhang et al. (2005) on their 
deep, high resolution spectrum of the very high excitation planetary 
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nebula NGC 7027. This spectrum covers the entire range between 3300 
and 9130 A. 

Table 1.1 gives the electronic configurations and ionization potentials 
of the C, N, O, Ne, S, CI and Ar ions. Those ions corresponding to 
the same isoelectronic sequence have been listed in the same column. 
Such ions have similar families of lines and provide the same plasma 
diagnostics. Note that Hll regions contain only traces (if any) of ions 
with ionization potential larger than 54 eV. On the other hand, plane- 
tary nebulae, which are generally ionized by hotter stars, can have an 
important fraction of Ne v and ions of similar ionization potentials. 
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Table 1.1. Electronic configurations and ionization potentials (in eV) 
of ions of C, N, O, Ne, S, CI, Ar with ionization potentials up to 

126 eV. 

Tables 1.2-1.8 list, for the ions of C, N, O, Ne, S, CI, Ar respectively, 
the vacuum wavelength of each line (in A), the identification, the type 
of transition (TT) , the electronic configurations for the lower and upper 
levels, the spectroscopic terms for the lower and upper level (J- J), the 
transition probability, and the energies of the lower and upper level. 

Note that it is more common to use the wavelength in the air than the 
wavelength in the vacuum. The IAU standard for conversion from air to 
vacuum wavelengths is given in Morton (1991). For vacuum wavelengths 
A a in Angstroms, convert to air wavelength A a via: A a = A v /(1.0 + 
2.735182 x 10~ 4 + 131.4182/A? + 2.76249 x 10 8 /A*). 
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The identification in the tables is represented as: Fe I for allowed 
transitions, Fe I] for intcrcombination transitions, [Fe I] for forbidden 
transitions. 

The symbols for the transition types are: El: allowed transitions; Ml: 
magnetic dipole forbidden transitions; E2: electric quadrupole forbidden 
transitions. 

The transition probabilities have been taken from the Atomic Line List 
of Peter van Hoof, or from the Chianti Atomic Database at 
|http://www.arcetri.astro.it/science/chianti/chianti.html| if followed by 
an asterisk. 



A vac (A) Ion TT Terms J- J A ki (s _1 ) Levels (cm -1 ) 



4622.864 


[Ci] 


Ml 


3p 


- x s 


1 - 







16.40 - 


21648.01 


8729.52 


[Ci] 


E2 


l D 


- x s 


2 - 







10192.63 - 


21648.01 


9811.01 


[Ci] 


E2 


3p 


- r D 


- 


2 




0.00 - 


10192.63 


9826.82 


[Ci] 


Ml 


3p 


- X D 


1 - 


2 




16.40 - 


10192.63 


9852.96 


[Ci] 


Ml 


3p 


- X D 


2 - 


2 




43.40 - 


10192.63 


2304000. 


[Ci] 


E2 


3p 


- 3 P 


- 


2 




0.00 - 


43.40 


3704000. 


[Ci] 


Ml 


3p 


- 3 P 


1 - 


2 


2.65(-7) 


16.40 - 


43.40 


6100000. 


[Ci] 


Ml 


3p 


- 3 P 


- 


1 


7.93(-8) 


0.00 - 


16.40 


1334.5323 


Cn 


El 


2po 


- 2 D 


i 

2 


3 
2 


2.38(+8) 


0.00 - 


74932.62 


1335.6627 


Cn 


El 


2po 


- 2 D 


3 
2 


3 
2 


4.74(+7) 


63.42 - 


74932.62 


1335.7077 


Cn 


El 


2po 


- 2 D 


3 
2 


o 
« 
2 


2.84(+8) 


63.42 - 


74930.10 


2322.69 


[Cn] 


M2 


2po 


- 4 P 


1 

i " 
2 




0.00 - 


43053.60 


2324.21 


Cn] 


El 


2po 


_4p 


2 




0.00 - 


43025.30 


2325.40 


Cn] 


El 


2po 


- 4 P 


1 

2 


1 
2 




0.00 - 


43003.30 


2326.11 


Cn] 


El 


2po 


_4p 


3 
2 




63.42 - 


43053.60 


2327.64 


Cn] 


El 


2po 


_4p 


3 


3 

i 




63.42 - 


43025.30 


2328.84 


Cn] 


El 


2po 


_4p 


3 
i 

i ' 

2 




63.42 - 


43003.30 


1576800. 


[Cn] 


Ml 


2pO 


_ 2po 




2.290(-6)* 


0.00 - 


63.42 


4550000. 


[Cn] 


Ml 


*P 


_4p 


2 


2.39(-7) 


43003.30 - 


43025.30 


1906.683 


[Cm] 


M2 




_ 3po 


- 


2 




0.00 - 


52447.11 


1908.734 


Cm] 


El 


*s 


_ 3po 


- 


1 




0.00 - 


52390.75 


1249200. 


[Cm] 


E2 


3po 


_ 3po 


- 


2 




52367.06 - 


52447.11 


1774300. 


[Cm] 


Ml 


3po 


_ 3po 


1 - 


2 


2.450(-6)* 


52390.75 - 


52447.11 


4221000. 


[Cm] 


Ml 


3po 


_ 3po 


- 


1 


2.39(-7) 


52367.06 - 


52390.75 



1548.203 Civ El 2 S - 2 P° |-| 2.65(+8) 0.00 - 64591.00 

1550.777 Civ El 2 S - 2 P° |-| 2.64(+8) 0.00 - 64483.80 



Table 1.2. Lines of C ions. 
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Avbc (A) Ion TT Terms J- J A ki (s 1 ) Levels (cm x ) 



3467.4898 


[Ni] 


Ml 


4 S° 


_ 2po 


3 
1 


3 

? 

2 


6.210(- 


-3)* 


0.00 - 


28839.31 


3467.5362 


[Ni] 


Ml 


4 S° 


_ 2po 


2.520(- 


-3)* 


0.00 - 


28838.92 


5199.3490 


[Ni] 


E2 


4 S° 


- 2 D° 


'i 
2 


3 
2 


2.260(- 


-5)* 


0.00 - 


19233.18 


5201.7055 


[Ni] 


E2 


4 S° 


- 2 D° 


3 


5 

§ 

2 


5.765(- 


-6)* 


0.00 - 


19224.46 


10400.587 


[Ni] 


Ml 


2 D° 


_ 2po 


2 


5.388(- 


-2)* 


19224.46 - 


28839.31 


10401.004 


[Ni] 


E2 


2 D° 


_ 2po 


2 


1 
2 


3.030(- 


-2)* 


19224.46 - 


28838.92 


10410.021 


[Ni] 


Ml 


2 D o 


_ 2po 


3 
2 


3 
2 


2.435(- 


-2)* 


19233.18 - 


28839.31 


10410.439 


[Ni] 


Ml 


2 D° 


_ 2po 


3 
2 


1 

2 


4.629(- 


-2)* 


19233.18 - 


28838.92 





r n tt i 


M2 


3p _ 


5 gO 


{} 


_ 2 






00 — 


46784 56 


91 39 683 


N II ] 


El 


3p _ 


5 go 


1 


_ 2 






48 67 — 


46784 56 


oi /l ^ /l 

Z 1410 .^OU 


N II ] 


E 1 


3p _ 


5 Q o 
O 


9 


_ 2 






i "^n sn 


/lfi7S/l 


3063.728 


[Nil] 


Ml 


3p_ 


\s 


1 


- 


4.588( 


-6)* 


48.67 - 


32688.64 


3071.457 


[Nil] 


E2 


3p_ 




2 


- 






130.80 - 


32688.64 


5756.240 


[Nil] 


E2 


1 D- 




2 


- 


9.234( 


-1)* 


15316.19 - 


32688.64 


6529.04 


[Nil] 


E2 


3p_ 







- 2 


1.928(- 


-6)* 


0.00 - 


15316.19 


6549.85 


[Nil] 


Ml 


3p_ 


X D 


1 


- 2 


9.819( 


-4)* 


48.67 - 


15316.19 


6585.28 


[Nil] 


Ml 


3p_ 


X D 


2 


- 2 


3.015( 


-3)* 


130.80 - 


15316.19 


764500. 


[Nil] 


E2 


3p_ 


3p 





- 2 






0.00 - 


130.80 


1217600. 


[Nil] 


Ml 


3p_ 


3p 


1 


- 2 


7.460(- 


-6)* 


48.67 - 


130.80 


2055000. 


[Nil] 


Ml 


3p_ 


3p 





- 1 


2.080( 


-6)* 


0.00 - 


48.67 


1744.351 


[Nm] 


M2 


2po _ 


4p 


i 

2 


5 
2 






0.00 - 


57327.90 


1746.823 


Nm] 


El 


2po _ 


4p 


1 
2 


3 
2 






0.00 - 


57246.80 


1748.646 


Nm] 


El 


2po _ 


4p 


1 

2 


1 
2 






0.00 - 


57187.10 


1749.674 


Nm] 


El 


2po _ 


4p 


3 
2 


5 
2 






174.40 - 


57327.90 


1752.160 


Nm] 


El 


2po _ 


4p 


3 
2 


3 

? 
2 






174.40 - 


57246.80 


1753.995 


Nm] 


El 


2po 


4p 


3 

i 

2 






174.40 - 


57187.10 


573400. 


[Nm] 


Ml 


2po _ 


2po 


3 

1 


4.770(- 


-5)* 


0.00 - 


174.40 


710000. 


[Nm] 


E2 


4p_ 


4p 


1 






57187.10 - 


57327.90 


1233000. 


[Nm] 


Ml 


4p _ 


4p 


'i 
I 

2 


1 


8.63(- 


6) 


57246.80 - 


57327.90 


1675000. 


[Nm] 


Ml 


4p _ 


4p 


§ 

2 


4.78(- 


6) 


57187.10 - 


57246.80 


1483.321 


[Niv] 


M2 




3pO 





- 2 






0.00 - 


67416.30 


1486.496 


Niv] 


El 




3pO 





- 1 






0.00 - 


67272.30 


482900. 


[Niv] 


E2 


3po 


3pO 





- 2 






67209.20 - 


67416.30 


694000. 


[Niv] 


Ml 


3po _ 


3po 


1 


- 2 


4.027( 


-5)* 


67272.30 - 


67416.30 


1585000. 


N iv 


Ml 


3po _ 


3po 





- 1 


4.52(- 


6) 


67209.20 - 


67272.30 



Table 1.3. Lines of N ions. 
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^vac (A) 


Ion 


TT 


Terms 


J - J 


A ki (s L ) 


Levels ( 


cm ) 


2973.1538 


[Oi] 


Ml 


3p _ 




1-0 




158.26 - 


33792.58 


5578.8874 


[Oi] 


E2 


1 D - 


lg 


2-0 




15867.86 - 


33792.58 


6302.046 


[Oi] 


Ml 


3p_ 


X D 


2-2 




0.00 - 


15867.86 


6365.536 


[Oi] 


Ml 


3p _ 


1 D 


1-2 




158.26 - 


15867.86 


6393.500 


[Oi] 


E2 


3p _ 


X D 


0-2 




226.98 - 


15867.86 


440573. 


[Oi] 


E2 


3p _ 


3p 


2-0 




0.00 - 


226.98 


631850. 


[Oi] 


Ml 


3p _ 


3p 


2-1 


8.91(-5) 


0.00 - 


158.26 


1455350. 


[Oi] 


Ml 


3p _ 


3p 


1-0 


1.75(— 5) 


158.26 - 


226.98 


2470.966 


[On] 


Ml 


4 s o _ 


2po 


f ~ I 


2.380(-2)* 


0.00 - 


40470.00 


2471.088 


[On] 


Ml 


4 S° - 


2po 




5.800(-2)* 


0.00 - 


40468.01 


3727.092 


[On] 


E2 


4 s o _ 


2 D° 


I I 


1.810(-4)* 


0.00 - 


26830.57 


3729.875 


[On] 


E2 


4 S° - 


2 D° 




3.588(-5)* 


0.00 - 


26810.55 


7320.94 


[On] 


E2 


2 D° - 


2po 


I " f 


5.630(-2)* 


26810.55 - 


40470.00 


7322.01 


[On] 


Ml 


2 D° - 


2po 


I I 


1.074(-1)* 


26810.55 - 


40468.01 


7331.68 


[On] 


Ml 


2 D o _ 


2po 


1 1 


9.410(-2)* 


26830.57 - 


40470.00 


7332.75 


[On] 


Ml 


2 D° - 


2po 


1 " 1 


5.800(-2)* 


26830.57 - 


40468.01 


4995000. 


[On] 


Ml 


2 D° - 


2 D° 


1 I 

2 2 


1.320(-7)* 


26810.55 - 


26830.57 


1657.6933 


[Om] 


M2 


3p _ 


5 S° 


0-2 




0.00 - 


60324.79 


1660.8092 


Ora] 


El 


3p _ 


5go 


1-2 




113.18 - 


60324.79 


1666.1497 


Om] 


El 


3p _ 


5 S° 


2-2 




306.17 - 


60324.79 


2321.664 


[Om] 


Ml 


3p _ 


x s 


1-0 




113.18 - 


43185.74 


2332.113 


[Om] 


E2 


3p _ 


x s 


2-0 




306.17 - 


43185.74 


4364.436 


[Om] 


E2 


X D - 


x s 


2-0 


4.760(-l)* 


20273.27 - 


43185.74 


4932.603 


[Om] 


E2 


3p _ 


1 D 


0-2 


7.250(-6)* 


0.00 - 


20273.27 


4960.295 


[Om] 


Ml 


3p_ 


X D 


1-2 


6.791(-3)* 


113.18 - 


20273.27 


5008.240 


[Oml 

L w J 


Ml 


3p _ 


X D 


2-2 


2.046(-2)* 


306.17 - 


20273.27 


326612. 


[Om] 


E2 


3p _ 


3p 


0-2 




0.00 - 


306.17 


518145. 


[Om] 


Ml 


3p _ 


3p 


1-2 


1.010( — 4)* 


113.18 - 


306.17 


883560. 


[Om] 


Ml 


3p_ 


3p 


0-1 


2.700( — 5)* 


0.00 - 


113.18 


1393.621 


[Oiv] 


M2 


2po 


4p 


1 5 




0.00 - 


71755.50 


1397.232 


Oiv] 


El 


2po _ 


4p 


? I 




0.00 - 


71570.10 


1399.780 


Oiv] 


El 


2po _ 


4p 


1 f 

2 2 




0.00 - 


71439.80 


1401.164 


Oiv] 


El 


2po _ 


4p 


3 5 
2 2 




386.25 - 


71755.50 


1404.813 


Oiv] 


El 


2po _ 


4 P 


3 3 




386.25 - 


71570.10 


1407.389 


Oiv] 


El 


2po _ 


4p 


I I 




386.25 - 


71439.80 


258903. 


[Oiv] 


Ml 


2po 


2po 


1 I 

2 2 


5.310(-4)* 


0.00 - 


386.25 


316800. 


[Oiv] 


E2 


4p_ 


4p 


1 5 

2 3 




71439.80 - 


71755.50 


539400. 


[Oiv] 


Ml 


4p_ 


4p 


3 5 
2 2 


1.03(-4) 


71570.10 - 


71755.50 


767000. 


[Oiv] 


Ml 


4p_ 


4p 


1 3 

2 2 


4.97(-5) 


71439.80 - 


71570.10 


1213.809 


[Ov] 


M2 


x s- 


3po 


0-2 




0.00 - 


82385.30 


225800. 


[Ov] 


E2 


3po _ 


3po 


0-2 




81942.50 - 


82385.30 


326100. 


[Ov] 


Ml 


3po 


3po 


1-2 


3.913(-4)* 


82078.60 - 


82385.30 


735000. 


[Ov] 


Ml 


3po _ 


3po 


- 1 


4.53(-5) 


81942.50 - 


82078.60 



Table 1.4. Lines of O ions. 
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\ ( X 

^vac {-^} 


Ion 


TT 


Terms 


J - 


J 


Aki (s 


) 


Levels ( 


„, — l \ 

cm ) 


128135.48 


[Nell] 


Ml 


2po 


2po 


3 

2 _ 


i 

2 


8.590(- 


3)* 


0.00 - 


780.42 


1793.802 


[ Ne in ] 


E2 


3p _ 


x s 


2 - 









0.00 - 


55747.50 


1814.730 


[ Ne in ] 


Ml 


3p _ 


x s 


1 - 









642.88 - 


55747.50 


3343.50 


[ Ne III ] 


E2 


1 D - 


x s 


2 - 





2.550 




25838.70 - 


55747.50 


3870.16 


[ Ne III ] 


Ml 


3p _ 


X D 


2 - 


2 


1.653(- 


1)* 


0.00 - 


25838.70 


3968.91 


[ Ne III ] 


Ml 


3p _ 


1 D 


1 - 


2 


5.107(- 


2)* 


642.88 - 


25838.70 


4013.14 


[ Ne III ] 


E2 


3p _ 


1 D 


- 


2 


1.555(- 


5)* 


920.55 - 


25838.70 


108630.5 


[ Ne III ] 


E2 


3p _ 


3p 


2 - 





2.524(- 


8)* 


0.00 - 


920.55 


155550.5 


[ Ne III ] 


Ml 


3p _ 


3p 


2 - 


1 


5.971(- 


3)* 


0.00 - 


642.88 


360135. 


[ Ne III ] 


Ml 


3p_ 


3p 


1 - 





1.154(- 


3)* 


642.88 - 


920.55 


4713.2 


[ Ne IV ] 


Ml 


2 D o _ 


2po 


5 
2 


3 
2 


4.857(- 


1)* 


41210. - 


62427. 


4715.4 


[ Ne iv j 


E2 


2 D o _ 


2po 


5 
2 


1 
2 


1.142(- 


1)* 


41210. - 


62417. 


4718.8 


[ Ne iv j 


Ml 


2 D o _ 


2po 


3 
2 


3 


5.881(- 


1)* 


41235. - 


62427. 


4721.0 


[ Ne iv j 


Ml 


2 D o _ 


2po 


3 
2 


I 
2 


4.854(- 


1)* 


41235. - 


62417. 


1575.1 


[ Ne v ] 


Ml 


3p _ 


x s 


1 - 









411.23 - 


63900.00 


1592.6 


[Nev] 


E2 


3p _ 


x s 


2 - 









1109.47 - 


63900.00 


2976. 


[Nov] 


E2 


r D - 


x s 


2 - 





2.834 




30294.00 - 


63900.00 


3301.0 


[Nov] 


E2 


3p_ 


X D 


- 


2 


5.460(- 


5)* 


0.00 - 


30294.00 


3346.4 


[Nev] 


Ml 


3p_ 


X D 


1 - 


2 


1.222(- 


1)* 


411.23 - 


30294.00 


3426.5 


[Nov] 


Ml 


3p_ 


X D 


2 - 


2 


3.504(- 


1)* 


1109.47 - 


30294.00 


90133.3 


[Nov] 


E2 


3p_ 


3p 


- 


2 






0.00 - 


1109.47 


143216.8 


[Nov] 


Ml 


3p_ 


3p 


1 - 


2 


4.585(- 


3)* 


411.23 - 


1109.47 


243175. 


[Nov] 


Ml 


3p _ 


3p 


- 


1 


1.666(- 


3)* 


0.00 - 


411.23 



Table 1.5. Lines of Ne ions. 
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What can emission lines tell us? 



^vac (A) 


Ion 


TT 


Terms 


j - j 




Levels (cm ) 


4590.5464 


[Si] 


Ml 


3 P - 1 S 


1-0 




396.06 - 22179.95 


7727.172 




1 

1 


E2 


X D - X S 


2-0 




9238.61 - 22179.95 


10824.140 


r o 


1 

J 


Ml 


3 P - 1 D 


2-2 




0.00 - 9238.61 


11308.950 


[S 


] 


Ml 


3 P - 1 D 


1-2 




396.06 - 9238.61 


11540.722 


[S 


i 


E2 


3 P - 1 D 


0-2 




573.64 - 9238.61 


174325.4 


r ^ 
L G 


1 
J 


E2 


l P ~ l P 


2-0 




0.00 - 573.64 


252490. 


r <s 
L G 


1 
J 


Ml 




2-1 


1.40(-3) 


0.00 - 396.06 


563111. 




1 
J 


Ml 


3p _ 3p 


1-0 


3.02(-4) 


396.06 - 573.64 


1250.5845 


s 


II 


El 




f ~ I 


4.43(+7) 


0.00 - 79962.61 


1253.8111 


c 
o 


II 


El 


S° - P 


1 I 

2 2 


4.40(+7) 


0.00 - 79756.83 


1259.5190 


Q 

o 


II 


El 


4 S o _ 2 P o 


4.34(+7) 


0.00 - 79395.39 


4069.749 


r q 


J 1 


Ml 




I _ I 


2.670(-l)* 


0.00 - 24571.54 


4077.500 


r q 




Ml 


S° - P° 


1.076(-1)* 


0.00 - 24524.83 


6718.29 


r a 




E2 


4 S° - 2 D° 


1 I 
1 I 
1 I 

1 I 

2 2 


3.338(-4)* 


0.00 - 14884.73 


6732.67 


r q 




E2 


i S° - 2 D° 


1.231(-3)* 


0.00 - 14852.94 


10289.55 


[S 


I] 


Ml 


o 2 p o 


1.644(-1)* 


14852.94 - 24571.54 


10323.32 


[S 


I] 


Ml 


2t^O 2nO 

Z D - P 


1.938(-1)* 


14884.73 - 24571.54 


10339.24 


r <=i 
L G 


.1 
1 J 


Ml 


2t^O 2nO 

Z D - P 


2 — 2 


1.812(-1)* 


14852.94 - 24524.83 


10373.34 


r 9 
L a 


.1 
1 J 


E2 


2 o 2 p o 


2^2 


7.506(-2)* 


14884.73 - 24524.83 


2141000. 


r c; 
L a 


1 J 


Ml 


2po _ 2po 


1 3 

2 2 


9 16C— 71 


24524.83 - 24571.54 


3146000. 


r c; 
L a 


1 J 


Ml 


2 D o _ 2 D v 


3 5 
2 2 


3.452(-7)* 


14852.94 - 14884.73 


1704.3928 


r a 


III 1 

III] 


M2 


3 J3 5 g o 


0-2 




0.00 - 58671.92 


1713.1137 


Q 

o 


TTT 1 

III] 


El 


3 j> 5 g o 


1-2 




298.68 - 58671.92 


1728.9415 


Q 

a 


TTI 1 

III J 


El 


3 j> 5 g o 


2-2 




833.06 - 58671.92 


3722.69 


r q 


TT 1 


Ml 


3p _ l g 


1-0 


1.016 


298.68 - 27161.00 


3798.25 


[S 


III] 


E2 


3 p _ lg 


2-0 




833.06 - 27161.00 


6313.8 


[S 


II ] 


E2 


D - S 


2-0 


3.045 


11322.70 - 27161.00 


8831.8 




TT 1 


E2 


3 P - 1 D 


0-2 


2.570(~4)* 


0.00 - 11322.70 


9071.1 


r 9 
L a 




Ml 


3 P - 1 D 


1-2 


3.107(-2)* 


298.68 - 11322.70 


9533.2 




TT 1 


Ml 


3 P - X D 


2-2 


1.856(-1)* 


833.06 - 11322.70 


120038.8 


[S 


II] 


E2 


3p _ 3p 


0-2 


5.008(-8)* 


0.00 - 833.06 


187130.3 


[S 


II] 


Ml 


P ~ P 


1-2 


1.877(— 3)* 


298.68 - 833.06 


334810. 


[S 


II] 


Ml 


3 T-i 3n 
P ~ P 


0-1 


8.429(-4)* 


0.00 - 298.68 


1387.459 


[S 


IV] 


M2 


2po _ 4p 


1 5 




0.00 - 72074.20 


1398.050 


s 


IV] 


El 


2po _ 4p 


1 § 
? i 

2 2 




0.00 - 71528.20 


1404.800 


s 


IV] 


El 


2pO _ 4p 




0.00 - 71184.50 


1406.019 


s 


IV] 


El 


2po _ 4p 


3 5 




951.43 - 72074.20 


1416.897 


s 


IV] 


El 


2pO _ 4p 


1 I 

2 2 




951.43 - 71528.20 


1423.831 


s 


IV] 


El 


2po _ 4p 


3 1 

? § 




951.43 - 71184.50 


105104.9 


[S 


v] 


Ml 


2po _ 2po 


7.750(-3)* 


0.00 - 951.43 


112400. 


[S 


IV] 


E2 


4p _ 4p 


? 1 

1 1 

2 2 




71184.50 - 72074.20 


183200. 


[S 


IV] 


Ml 


4p _ 4p 


2.63(-3) 


71528.20 - 72074.20 


291000. 


[S 


IV] 


Ml 


4p_4p 


1 3 

2 2 


9.13(-4) 


71184.50 - 71528.20 


1188.281 


[S 


v] 


M2 


lg _ 3 p o 


0-2 




0.00 - 84155.20 


88400. 


[S 


v] 


E2 


3pO _ 3pO 


0-2 




83024.00 - 84155.20 


131290. 


[S 


v] 


Ml 


3po _ 3po 


1-2 


5.96(-3) 


83393.50 - 84155.20 


270600. 


[S 


v] 


Ml 


3po _ 3po 


- 1 


9.07(-4) 


83024.00 - 83393.50 



Table 1.6. Lines of S ions. 
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A vac 


Ion 


TT 


Terms 


J 


— J 


Aki (s ) 


Levels ( 


cm ) 


113333.52 


[Cli] 


Ml 


2po 


- 2 P° 


3 
2 


1 

~ 2 


1.24( — 2) 


0.00 - 


882.35 


3587.055 


[CI ii] 


E2 


3p 


- x s 


2 


- 




0.00 - 


27878.02 


3678.902 


[CI ii] 


Ml 


3p 


- x s 


1 


- 


1.297 


696.00 - 


27878.02 


6163.54 


[CI ii] 


E2 


X D 


- l s 


2 


- 


2.416 


11653.58 - 


27878.02 


8581.05 


[CI ii] 


Ml 


3p 


- 1 D 


2 


- 2 


1.038(-1)* 


0.00 - 


11653.58 


9126.10 


[CI ii] 


Ml 


3p 


- l D 


1 


- 2 


2.872(-2)* 


696.00 - 


11653.58 


9383.41 


[CI ii] 


E2 


3p 


- X D 





- 2 


1.134(-5)* 


996.47 - 


11653.58 


100354. 


[CI ii] 


E2 


3p 


- 3 P 


2 


- 


6.053(-7)* 


0.00 - 


996.47 


143678. 


[CI ii] 


Ml 


3p 


- 3 P 


2 


- 1 


7.566(^3)* 


0.00 - 


696.00 


332810. 


[CI ii] 


Ml 


3p 


- 3 P 


1 


- 


1.462( — 3)* 


696.00 - 


996.47 


3343.81 


[CI in] 


Ml 


4 S° 


_ 2po 


§ 


~ § 




0.00 - 


29906.00 


3354.17 


[CI in] 


Ml 




_ 2po 




~ 2 




0.00 - 


29813.60 


5519.25 


[CI in] 


E2 


4 S° 


- 2 D° 


2 


~ | 




0.00 - 


18118.40 


5539.43 


[CI in] 


E2 


4 S° 


- 2 D° 


| 






0.00 - 


18052.40 


8436.3 


[CI in] 


Ml 




_ 2po 


2 


2 




18052.40 - 


29906.00 


8483.5 


[CI m] 


Ml 


2 D° 


_ 2po 


5 


3 




18118.40 - 


29906.00 


8502.5 


[CI in] 


Ml 


2 D° 


- 2 P° 


2 


? 
2 




18052.40 - 


29813.60 


8550.5 


[ CI in ] 


E2 


2 D° 




2 


i 

~ 2 




18118.40 - 


29813.60 


1082000. 


[ CI in ] 


Ml 


2po 


_ 2po 


1 

I 




7.09(-6) 


29813.60 - 


29906.00 


1515000. 


[CI in] 


Ml 


2 D o 


- 2 D° 


2 


" f 

_ 2 


3.10(-6) 


18052.40 - 


18118.40 


3119.51 


[CI iv] 


Ml 


3p 


- *S 


1 


- 




492.50 - 


32548.80 


3204.52 


[ CI iv ] 


E2 


3p 


- 1 S 


2 


- 




1342.90 - 


32548.80 


5324.47 


[ CI iv ] 


E2 


X D 


-'s 


2 


- 




13767.60 - 


32548.80 


7263.4 


[CI iv] 


E2 


3p 


- r D 





- 2 




0.00 - 


13767.60 


7532.9 


[CI iv] 


Ml 


3p 


- r D 


1 


- 2 




492.50 - 


13767.60 


8048.5 


[CI iv] 


Ml 


3p 


- X D 


2 


- 2 




1342.90 - 


13767.60 


74470. 


[CI iv] 


E2 


3p 


_3p 





- 2 




0.00 - 


1342.90 


117590. 


[CI iv] 


Ml 


3p 


_3p 


1 


- 2 


8.29(-3) 


492.50 - 


1342.90 


203000. 


[CI iv] 


Ml 


3p 


_3p 





- 1 


2.15(-3) 


0.00 - 


492.50 


67090. 


[Civ] 


Ml 


2po 


_ 2po 


l 

2 


3 
2 


2.98(-2) 


0.00 - 


1490.50 


72400. 


[Civ] 


E2 


4p 


_4p 


1 
2 


~ 2 




86000. +v - 


87381. +v 


118600. 


[Civ] 


Ml 


4p 


_4p 


3 
2 


' f 


9.70(-3) 


86538. +v - 


87381. +v 


186000. 


[Civ] 


Ml 


4p 


_4p 


1 

2 


2 


3.50(-3) 


86000. +v - 


86538. +v 


58200. 


[ CI vi ] 


E2 


3po 


_ 3po 





- 2 




97405. + x - 


99123. + x 


85800. 


[ CI vi j 


Ml 


3po 


_ 3po 


1 


- 2 


2.13(-2) 


97958. + x - 


99123. + x 


180800. 


[ CI vi ] 


Ml 


3po 


_ 3po 





- 1 


3.04(-3) 


97405. + x - 


97958. + x 



Tabic 1.7. Lines of CI ions. 
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\ ( X 

^vac {-^} 


Ion 


TT 


Terms 


J - 


J 


A-ki (s ) 


Levels ( 


cm ) 


69852.74 


[Aril] 


Ml 


P° - P° 


3 

2 _ 


1 

2 


5.28(-2) 


0.00 - 


1431.58 


3006.10 


[ Ar III ] 


E2 


3p _ lg 


2 - 







0.00 - 


33265.70 


3110.08 


[ Ar III j 


Ml 


3 P - l s 


1 - 







1112.18 - 


33265.70 


5193.27 


[ Ar III j 


E2 


1 D - 1 S 


2 - 







14010.00 - 


33265.70 


7137.8 


[ Ar III ] 


Ml 


3 P - 1 D 


2 - 


2 




0.00 - 


14010.00 


7753.2 


[Arm] 


Ml 


3 P - 1 D 


1 - 


2 




1112.18 - 


14010.00 


8038.8 


[ Ar in ] 


E2 


3 P - *D 


- 


2 




1570.28 - 


14010.00 


63682.9 


[Arm] 


E2 


?>p _ 3p 


2 - 







0.00 - 


1570.28 


89913.8 


[ Ar in ] 


Ml 


3 f ~ 3 P 


2 - 


1 


3.09(-2) 


0.00 - 


1112.18 


218291. 


[ Ar in ] 


Ml 


P - P 


1 - 





5.19(-3) 


1112.18 - 


1570.28 


2854.48 


[ Ar iv ] 


Ml 


4 S° - 2 P° 


| ~ 




2.134 


0.00 - 


35032.60 


2868.99 


j Ar IV j 


Ml 


S° - P° 


2 ~ 


Y 

2 


8.702(-l)* 


0.00 - 


34855.50 


4712.69 


[ Anv] 


E2 


S° - D° 


2 ~ 




1.906(-3)* 


0.00 - 


21219.30 


4741.49 


[ Ar IV ] 


E2 


A S° - 2 D° 


§ " 


i 


2.277(-2)* 


0.00 - 


21090.40 


7172.5 


[ Ar IV ] 


Ml 


9 n 9 
^T) J ^JD iJ 


2 


1 

2 


8.915(-1)* 


21090.40 - 


35032.60 


7239.4 


[ Ar IV j 


Ml 


2jjO _ 2pO 


5 


3 


6.559(-l)* 


21219.30 - 


35032.60 


7264.7 


[Anv] 


Ml 


2 D° - 2 P° 


§ 
2 


f 
2 


6.660(-l)* 


21090.40 - 


34855.50 


7333.4 


[ Ar iv ] 


E2 


2 D° - 2 P° 


2 ~ 


; 

2 


1.179(-1)* 


21219.30 - 


34855.50 


564700. 


[ Ar iv ] 


Ml 


2p o 2 p o 


1 

I " 


3 
2 


4.937( — 5)* 


34855.50 - 


35032.60 


776000. 


[ Anv] 


Ml 


D - D 


2 


2 


2.287(-5)* 


21090.40 - 


21219.30 


2691.63 


[ Ar v 1 


Ml 


3p _ lg 


1 - 


o 




763.23 - 


37915.50 


2786.55 


j Ar v] 


E2 


3p_lg 


2 - 







2028.80 - 


37915.50 


4626.22 


[ Ar v ] 


E2 


D - S 


2 - 







16299.60 - 


37915.50 


6135.12 


[ Ar v] 


E2 


3 P - 1 D 


- 


2 




0.00 - 


16299.60 


6436.51 


[ Ar v] 


Ml 


3 P - 1 D 


1 — 


2 




763.23 - 


16299.60 


7007.3 


[ Ar v] 


Ml 


3p _ 1 D 


2 - 


2 




2028.80 - 


16299.60 


49290.2 


[ Ar v] 


E2 


3p _ 3 p 


- 


2 




0.00 - 


2028.80 


79015.8 


j Ar v] 


Ml 


3p _ 3p 


1 - 


2 


2.73(-2) 


763.23 - 


2028.80 


131021.9 


[ Ar v] 


Ml 


3p _ 3p 


- 


1 


7.99(-3) 


0.00 - 


763.23 


978.58 


[ Ar vi ] 


M2 


2po _ 4p 


1 

2 


5 
2 




0.00 - 


102189.20 


45295. 


j Ar vi j 


Ml 


2po _ 2po 


1 

2 


3 
2 


9.68(-2) 


0.00 - 


2207.74 


49237. 


[Arvij 


E2 


4p_4p 


1 

2 


5 




100158.20 - 


102189.20 


81010. 


[ Ar VI ] 


Ml 


4p _ 4p 


3 
2 


2 


3.04(-2) 


100954.80 - 


102189.20 


125530. 


[ Ar vi ] 


Ml 


4p_4p 


1 

2 


3 
2 


1.14(-2) 


100158.20 - 


100954.80 


865.16 


[ Ar VII ] 


M2 


lg_3po 


- 


2 




0.00 - 


115585.00 


40310. 


j Ar VII ] 


E2 


3 p o _ 3 p o 


- 


2 




113104.00 - 


115585.00 


59500. 


[ Ar VII ] 


Ml 


3 p o _ 3 p o 


1 - 


2 


6.41(-2) 


113904.00 - 


115585.00 


125000. 


j Ar VII ] 


Ml 


3 p o _ 3 p o 


- 


1 


9.21(-3) 


113104.00 - 


113904.00 



Table 1.8. Lines of Ar ions. 



In Tables 1.2-1.8, it is easy to see which ions are available for each 
element in a given wavelength range. For example, there is no line of N in 
in the optical. One can also see which lines arise from the same upper 
level, and therefore have intensity ratios that depend only on atomic 
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physics and not on the conditions inside the nebula. For example, for 
[O in] , this is the case for the 4969 and 5007A lines, whose intensity ratio 
is given by the ratio of transition probabilities divided by the ratio of 
wavelengths. 

Tables 1.9-1.14 group the lines according to the electronic configu- 
ration of their ions. Fig. 1.6 shows simplified energy level diagrams, 
which allow one to understand which line pairs are temperature or den- 
sity indicators. Of course, in practice, one uses only the lines with largest 
transition probabilities. For example, [Oin] A4363/4931 could be a tem- 
perature indicator as [Oni] A4663/5007 but the intensity of the A4931 
line is less than one thousand's of that of A5007 (as inferred from Table 
1.4, comparing the transition probabilities). 



T Ai 


> 




lis 


2 |s 


3 








, R R R. 


V 


\ 






\ 




< * 



1 So 

1 D 2 



3 Po,1,2 



1 Ai A2A3A4A5R _, 

I=] 2 P 



N1N2 



± ± 



1/2,3/2 



^J 2 D 3 /2,5 



Fig. 1.6. Simplified energy level diagrams with transition types indicated. A stands 
for "auroral", N for "nebular", T for "transauroral" and F for "fine structure"". Left: 
ions p 4 (Table 1.10) and p 2 (Table 1.12). Right: ions p 3 (Table 1.11) 





[Ne 11] 


[Cli] 


[Aril] 


Fi 


128135.48 


113333.52 


69852.74 



Table 1.9. Lines from 2p 5 and 3p 5 ions. The first column gives the 
type of line, as defined in Fig. 1.6. 
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[Oi] 


[Nc in] 


[Si] 


[CI II] 


[ Ar III ] 






Ti 


2973.1538 


1814.730 


4590.5464 


3678.902 


3110.08 






Ai 


5578.8874 


3343.50 


7727.172 


6163.54 


5193.27 






Ni 


6302.046 


3870.16 


10824.140 


8581.05 


7137.8 






N 2 


6365.536 


3968.91 


11308.950 


9126.10 


7753.2 






N 3 


6393.500 


4013.14 


11540.722 


9383.41 


8038.8 






Fi 


440573. 


108630.5 


174325.4 


100354. 


63682.9 






F 2 


631850. 


155550.5 


252490. 


143678. 


89913.8 






F 3 


1455350. 


360135. 


563111. 


332810. 


218291. 






Table 1.10. Lines from 2p 4 


and 3p 4 ions. 


The first column gives the 








type of line, 


as defined in 


Fig. 1.6. 










[Ni] 


[On] 


[Noiv] 


;sn] 


[CI in] 


[ Ar IV ] 




Ti 


3467.4898 


2470.966 


4713.2 4069.749 


3343.81 


2854.48 




T 2 


3467.5362 


2471.088 


4715.4 4077.500 


3354.17 


2868.99 




Ni 


5199.3490 


3727.092 


4718.8 6718.29 


5519.25 


4712.69 




N 2 


5201.7055 


3729.875 


4721.0 6732.67 


5539.43 


4741.49 




Ai 


10400.587 


7320.94 


10289.55 


8436.3 


7172.5 




A 2 


10401.004 


7322.01 


10323.32 


8483.5 


7239.4 




A 3 


10410.021 


7331.68 


10339.24 


8502.5 


7264.7 




A 4 


10410.439 


7332.75 


10373.34 


8550.5 


7333.4 




Fi 




4995000. 


2141000. 1082000. 


564700. 




F 2 






3146000. 1515000. 


776000. 




Table 1.11. Lines from 2p 3 


and 3p 3 ions. 


The first column gives the 








type of line, 


as defined in 


Fig. 1.6. 










[Ci] 


[Nil] 


[Om] 


[Nov] 


[Sin] 


[CI iv] 


[Arv] 


Ti 


4622.864 


3063.728 


2321.664 


1575.1 


3722.69 


3119.51 


2691.63 


Ai 


8729.52 


5756.240 


4364.436 


2976. 


6313.8 


5324.47 


4626.22 


Ni 


9811.01 


6529.04 


4932.603 


3301.0 


8831.8 


7263.4 


6135.12 


N 2 


9826.82 


6549.85 


4960.295 


3346.4 


9071.1 


7532.9 


6436.51 


N 3 


9852.96 


6585.28 


5008.240 


3426.5 


9533.2 


8048.5 


7007.3 


Fi 


2304000. 


764500. 


326612. 


90133.3 


120038.8 


74470. 


49290.2 


F 2 


3704000. 


1217600. 


518145. 


143216.8 


187130.3 


117590. 


79015.8 


F 3 


6100000. 


2055000. 


883560. 


243175. 


334810. 


203000. 


131021.9 



Table 1.12. Lines from 2p 2 and 3p 2 ions. The first column gives the 
type of line, as defined in Fig. 1.6. 
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[Cn] 


[Nm] 


[Oiv] 


[Siv] 


[Civ] 


Fi 


1576800. 


573400. 


258903. 


105104.9 


67090. 


F 2 










72400. 


F 3 










118600. 


F 4 










186000. 



Table 1.13. Lines from 2p and 3p ions. The first column gives the 
type of line, as defined in Fig. 1.6. 





[Cm] 


[Niv] 


[Ov] 


Fi 


1249200. 


482900. 


225800. 


F 2 


1774300. 


694000. 


326100. 


F 3 


4221000. 


1585000. 


735000. 



Table 1.14. Lines from 2s 2 and 3s 2 ions. The first column gives the 
type of line, as defined in Fig. 1.6. 
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